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Abstract 


The  Department  of  Defense  has  identified  launch  vibration  isolation  as  a  major 
research  interest.  Reducing  the  loads  a  satellite  experiences  during  launch  will  greatly 
enhance  the  reliability,  lifetime,  and  payload  to  structure  ratio.  DoD  satellite  programs 
stand  to  benefit  significantly  from  advances  in  vibration  isolation  technology. 

This  study  explores  potential  hybrid  vibration  isolation  control  designs  versus 
passive  designs.  A  simple  lumped  mass  dynamic  model  of  a  satellite  and  a  representative 
launch  vehicle  was  designed  using  Simulink.  The  analysis  focuses  on  the  various  sources 
of  transient  launch  accelerations  such  as  aero-acoustic  loads,  separation  events,  wind 
gusts,  and  motor  induced  vibration.  The  passive  vibration  suppression  design  reduced 
axial  “bounce”  modes.  Further  reductions  were  made  possible  with  an  added  active 
controller. 

The  results  of  modeling  indicate  that  as  much  as  a  90  percent  improvement  in 
loads  on  the  satellite  were  recognized  from  the  combination  of  active  and  passive 
vibration  control  techniques.  The  model  also  explored  the  response  of  a  sub-component 
of  the  satellite  payload.  For  passive  isolation  the  sub-component  remained  near  its  low 
baseline  level  of  response,  but  for  hybrid  isolation  the  sub-component  is  exposed  to 
greater  levels  of  loading  as  a  function  of  passive  isolation  frequency. 
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Whole  Spacecraft  Vibration  Isolation 


1  Introduction 

In  almost  all  heavy  industry  endeavors,  vibration  suppression  or  isolation  is  a 
major  concern.  Preventing  earthquake  damage  by  mounting  high  rise  buildings  on 
massive  springs  and  quieting  the  ride  in  a  car  with  automotive  suspensions  are 
recognizable  examples  of  vibration  isolation  or  suppression.  Vibration  isolation  and 
suppression  allows  for  greater  comfort,  control  or  reduced  risk  depending  on  its 
application. 

Recently,  the  Air  Force  supported  the  Navy  GEOS  AT  Follow-on  satellite  launch 
on  the  Taurus  launch  vehicle.  The  Taurus  is  a  solid-propelled,  four-stage  launch  vehicle, 
which  is  capable  of  placing  small  satellites  into  low-earth  orbits.  The  satellite,  it  was 
shown  during  testing,  did  not  meet  the  required  standards  to  survive  the  launch 
environment  of  the  Taurus.  The  Air  Force  Research  Laboratory  (AFRL)  was  able  to 
provide  a  satellite  vibration  isolation  system  that  allowed  GFO  to  be  launched  as  is 
without  a  costly  redesign  [27]. 

This  thesis  seeks  to  extend  the  work  currently  being  undertaken  by  AFRL.  Using 
Simulink  and  Matlab  [3]  software  packages  several  models  were  developed  to  explore 
passive  and  hybrid  whole  satellite  vibration  isolation.  The  models  include  a  launch 
vehicle  dynamic  model  and  a  flight  model  to  re-create  the  loads  a  satellite  may  encounter 
during  a  launch.  The  research  explored  both  lateral  and  axial  vibration  modes  so  that  the 


1-1 


dynamic  coupling  of  the  isolation  system  can  be  analyzed.  Several  types  of  active  control 
systems  were  tried  in  the  model  to  demonstrate  the  different  benefits  and  shortcomings  of 
the  controllers.  The  final  model  used  was  designed  with  maximum  flexibility  so  that  the 
researcher  could  rapidly  change  system  parameters  and  explore  different  launch 
environmental  considerations. 

1.1  Background 

From  the  beginning  of  the  space  program  it  has  been  recognized  that  the  launch 
vehicle  vibration  environment  is  a  major  source  of  failures  of  satellites.  Over  time 
NASA,  the  DoD  and  commercial  launch  providers  have  developed  extensive  test 
requirements  for  satellites  in  an  effort  to  enhance  satellite  launch  survivability  from  the 
beginning  of  the  design  phase.  The  testing  program  is  costly  in  both  time  and  currency, 
and  frequently  leads  to  the  redesign  of  a  satellite. 

Personal  experience  has  shown  that  a  major  consideration  in  the  planning  stages 
of  a  launch  campaign  is  the  satellite  launch  environment  qualification  test  and 
verification  process.  This  process  includes  lengthy  discussions  of  the  satellite’s  design, 
testing  procedures,  and  coupled-loads  analysis  (CLA)  results.  Both  launch  provider  and 
launch  customer  often  require  CLA.  It  is  done  on  the  satellite  manufacturer’s  behalf  for 
assurances  that  the  design  is  sound  when  placed  in  the  environment  of  a  particular  launch 
vehicle  (LV).  The  process  is  repeated  by  the  launch  provider  so  that  the  dynamic 
response  of  the  satellite  to  a  LV  induced  loading  can  be  estimated  and  the  LV  guidance 
computer  programmed.  The  result  is  hundreds  of  hours  of  negotiation,  testing,  and 
modeling  to  confirm  the  safety  of  a  launch. 
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Early  references  to  satellite  failure  analysis  show  vibration  to  be  a  leading  source 
of  failure.  Goddard  Space  Flight  Center  conducted  a  study  of  first  day  launch  failures  in 
1971  [25].  This  was  folio  wed-up  three  years  later  with  an  analysis  of  first  month  failures 
[24],  Both  studies  are  of  the  same  57  spacecraft  launch  by  NASA  over  the  ten-year  period 
of  1960  to  1970.  The  one  LV  environmental  element  singled  out  in  the  conclusion  as  the 
source  of  30-60  percent  of  first  day  failures  is  the  vibration/acoustic  load.  The  paper  goes 
on  to  note  that  it  is  not  always  possible  to  identify  vibration  as  the  source  of  a  failure  once 
the  satellite  is  on  orbit,  but  that  ground  testing,  when  done  effectively,  very  often 
provides  evidence  that  the  source  of  a  particular  failure  is  vibration  related. 

The  Jet  Propulsion  Laboratory  has  developed  extensive  testing  procedures  to 
verify  the  validity  of  a  satellite  or  satellite  sub-component  design  prior  to  launch.  JPL 
also  dedicates  research  efforts  to  optimizing  future  spacecraft  testing.  Being  developed  is 
a  generic  launch  environment  envelope  that  would  allow  for  the  pre-qualification  of 
spacecraft  hardware,  thus  minimizing  testing  requirements.  This  testing  envelope 
attempts  to  encompass  95%  of  the  environmental  characteristics  of  nearly  all  current  US 
launch  vehicles  as  well  as  some  foreign  LVs  [12].  The  focus  of  recent  experimentation  at 
JPL  and  Los  Alamos  National  Laboratory  is  improving  the  realism  of  hardware  trials. 
Improving  realism  has  involved  creating  special  force  limiting  algorithms  to  notch  certain 
frequency  bands  near  resonant  frequencies  of  the  satellite  to  avoid  damaging  the  test 
article,  and  also  better  matching  the  characteristics  of  the  vibration  loads  anticipated 
during  launch.  The  key  difficulty  of  vibro-mechanical  testing  of  satellites  and 
components  is  that  the  shaker  has  infinite  impedance  and  the  satellite  clearly  does  not. 
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Matching  the  impedance  of  the  test  bench  to  the  test  article  improves  the  realism  of  a  test 
program  dramatically  [12],  [18],  [22]. 

Joshi  has  numerically  studied  the  effects  of  thrust-to-weight  ratio  versus 
transverse  stiffness.  The  work  looks  to  understand  the  relationship  of  trajectory  and  the 
fundamental  modes  of  a  launch  vehicle  structure.  The  key  finding  is  that  the  modes  of  a 
LV  are  constantly  changing  as  a  function  of  the  continuous  mass  change  of  the  rocket.  An 
overall  reduction  of  the  transverse  modes  results  due  to  the  increasing  axial  compression 
[7].  But  of  significance  to  this  study  is  the  constantly  changing  frequencies.  This 
variability  will  complicate  isolator  designs,  both  passive  and  hybrid. 

Another  environmental  factor  that  has  received  attention  is  wind  gust  buffet. 
Langley  Research  Center  has  conducted  wind-tunnel  tests  to  characterize  the  buffet 
response  of  hammerhead  launch  vehicles  [2].  The  research  sought  to  create  an  aero- 
elastically  scaled  model  of  the  upper  portion  of  the  Atlas  LV.  The  results  showed  there  is 
increased  buffet  due  to  the  configuration  of  the  payload  fairing,  but  the  ability  to 
accurately  characterize  the  buffet  mechanism  was  limited  by  the  wind-tunnel  test  facility. 
Prediction  of  environmental  loads  is  a  significant  part  of  launch  preparation.  The  launch 
vehicle  manufacturer  has  to  be  able  to  represent  the  environment  so  that  payload 
qualification  can  proceed  accurately. 

Coupled  loads  analysis  seeks  to  computationally  predict  the  response  of  a  payload 
to  the  flight  environment  at  various  stages  of  an  LVs  flight.  This  is  a  costly,  time 
consuming  process.  CLA  is  used  as  a  design  guide,  and  is  used  to  direct  testing  programs 
[23].  CLA  is  costly  because  of  the  time  required  to  create  and  verify  a  finite-element  or 


1-4 


statistical  energy  model  of  a  spacecraft  and  the  payload  interface.  The  results  may  be 
inaccurate  as  well. 

Efforts  by  Lee-Glauser,  Ahmadi  and  Layton  at  Langley  Research  Center  have 
focused  on  vibration  isolation  for  the  Space  Transportation  System  (STS).  They  have 
developed  passive  and  active  vibration  control  strategies  for  shuttle  payloads.  Using  a 
five  DoF  lump  parameter  model  the  satellite  payload  they  have  developed  several  isolator 
types.  The  models  are  used  to  simulate  the  isolation  effectiveness  of  the  first  eight 
seconds  of  flight.  For  the  STS,  this  is  the  period  of  greatest  vibration  energy  as  the  solid 
rocket  boosters  and  Space  Shuttle  Main  Engines  are  ignited  [13],  [14],  [15], 

As  mentioned  previously,  the  US  Air  Force  supported  the  launch  of  the  GEOSat 
Follow-on  Satellite  (GFO)  program  by  designing,  testing  and  delivering  a  passive  whole 
satellite  launch  vibration  isolation  system.  The  GFO  satellite  failed  to  qualify  for  launch 
in  its  pre-test  configuration  and  as  such  was  to  incur  a  launch  delay  to  improve  the  its 
response  to  test  loads  via  a  redesign.  This  delay  was  offset  by  the  design  of  a  simple 
passive  isolator  that  added  only  a  small  change  in  forward  location  of  the  satellite.  The 
Air  Force  contracted  the  design  out  to  CSA  Engineering.  CSA  rapidly  prototyped  the 
isolator  design  using  the  NASTRAN  FEM  software  package.  They  then  tested  the  design 
for  flight  qualification.  The  final  design  successfully  flew  with  GFO  on  February  10, 
1998  [27]. 

CSA  is  currently  designing  a  larger  payload  isolation  system  for  medium 
capability  launch  vehicles  and  3000-kg  class  satellites.  Their  design  constraints  for  a 
replacement  payload  attach  fitting  that  would  provide  lateral  vibration  isolation  are  as 
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follows.  The  isolator  could  not  create  structural  modes  below  6Hz  because  the  low 


frequency  modes  would  have  an  impact  on  the  LV  guidance,  control  and  navigation 
system.  The  second  constraint  was  to  limit  lateral  oscillations  to  less  than  one  inch  total 
displacement  or  violate  the  dynamic  fairing  envelope.  To  do  so  would  likely  cause  the 
failure  of  the  launch  if  the  satellite  struck  the  payload  fairing.  Their  work  on  the  GFO 
isolation  system  and  the  effort  described  here  is  part  of  a  broader  USAF  Air  Force 
Research  Lab  effort  in  whole  spacecraft  vibration  isolation  [28]. 

Air  Force  Research  Lab  granted  the  Small  Business  Innovative  Research  contracts 
to  CSA  for  their  work  on  whole  spacecraft  vibration  isolation.  The  Kirtland  AFB  efforts 
are  initial  experimental  work  on  hybrid  launch  isolation.  Their  initial  results  show  great 
gains  over  passive  isolation  systems  for  the  reduction  of  vibration  loads  of  a  test  article. 
Goals  they  have  determined  for  a  hybrid  isolator  are  fail-safe  design,  limited  weight 
increase,  adaptive  to  changing  environment,  broadband  and  narrowband  attenuation  of 
the  launch  dynamic  environment.  Achievement  of  these  goals  will  lead  to  decreased 
satellite  mass  (or  conversely  increased  useful  satellite  payload  mass),  increased  satellite 
reliability  and  survivability,  and  decreased  life-cycle  costs  [19]. 

As  a  final  note,  the  Department  of  Defense  has  indicated  its  top-level  interest  in 
vibration  isolation  studies  in  the  1998  Defense  Technology  Area  Plan  --  Space  Platforms 
DTOs.  The  goals  laid  out  by  DoD  are  the  reduction  of  launch  vehicle  structural  mass  by 
40%  and  satellite  launch  loads  80%  by  FY01.  These  two  are  conflicting  goals  because 
reduced  structural  mass  generally  results  in  higher  vibration  loading  on  the  satellite 


1-6 


payload.  Over  the  next  five  years  approximately  sixty  million  dollars  are  budgeted  for 
these  and  similar  efforts  to  be  performed  by  the  Air  Force  [1]. 

1.2  Research  Objectives 

Research  into  the  nature  of  the  dynamic  loads  on  a  satellite  has  been  pursued 
since  the  beginning  of  the  space  program.  Early  failure  analysis  points  to  vibro- 
mechanical  loads  as  the  key  failure  mode.  Characterization  of  the  loads  on  a  satellite  is 
on-going.  Passively  isolating  a  satellite  from  the  vibration  loads  created  by  the  launch 
environment  was  recently  achieved  with  the  GFO  satellite.  Efforts  to  further  control  the 
environmental  loads  on  a  satellite  are  currently  underway. 

While  CLA  can  explore  the  vibration  modes  of  a  satellite  and  payload  interface  it 
is  time-invariant,  focused  on  specific  load  cases  and  is  limited  in  its  application  to 
dynamic  control  design.  Though  highly  effective  as  a  satellite  qualification  tool,  CLA  is 
too  complex  for  the  analysis  of  the  dynamics  that  would  drive  a  hybrid  vibration  isolation 
system  design.  Therefore,  development  of  a  time-varying  model  of  a  launch  vehicle  and 
associated  environmental  loads  would  assist  the  design  process  of  a  hybrid  isolation 
system. 

This  study  discusses  the  researcher’s  development  of  a  time-variant  launch 
vehicle  and  launch  environment  model  upon  which  the  study  of  passive  and  hybrid 
vibration  isolation  methods  was  based.  The  researcher  discusses  the  results  of  analysis  on 
passive  and  hybrid  vibration  isolation  systems.  The  emphasis  of  the  isolation  research  is 
the  overall  reduction  of  whole  spacecraft  loads  and  a  payload  sub-component.  Finally, 
recommendations  for  future  efforts  in  this  area  are  summarized. 
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1.3  Organization  of  this  Paper 


This  paper  is  divided  into  three  sections.  The  first  is  background  information 
exploring  historical  and  current  efforts  in  launch  vibration  isolation.  The  second  section 
reviews  the  derivations  necessary  to  build  models  of  the  launch  vehicle  to  study  its 
vibration  characteristics,  discusses  the  passive  and  active  control  theories  pertinent  to  the 
design  of  passive  and  hybrid  isolation  systems,  and  then  overviews  the  top-level  design 
philosophy  of  the  Simulink  models  created  to  implement  the  equations  established  earlier 
in  the  section.  The  third  section  provides  the  analysis  of  passive  and  hybrid  simulation 
results.  The  appendices  will  provide  detailed  discussion  of  the  model  designs  and  source 
codes  of  the  analysis  tools  developed  for  this  research  effort. 
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2  Model  Development 


In  this  chapter  the  equations  used  to  describe  the  vibration  loads  and  other  forces 
that  act  on  the  launch  vehicle  and  its  payload  as  it  flies  to  orbit  are  developed.  The 
derivation  begins  with  the  fundamental  gravity  turn  equations.  The  axial  and  lateral 
vibration  equations  are  presented.  These  are  followed  by  brief  discussions  of  equations 
used  to  simulate  environmental  elements  such  as  vibration  loading.  The  last  section  of 
Chapter  2  details  the  integration  of  the  parts  into  a  Simulink  model  to  test  the  benefits  of 
passive  and  hybrid  isolation.  Detailed  discussion  of  the  simulation  sub-systems  is  in 
Appendix  A. 

2.1  Derivation  of  the  Equations  of  Motion 

Fundamental  to  the  generation  of  a  flight  environment  is  the  LY  trajectory  and  the 
quasi-static  loads  it  creates.  Quasi-static  loads  are  the  accelerations  generated  by  the 
relatively  slow  changes  of  thrust,  drag  and  gravity.  The  LV  trajectory  equations  are 
derived  in  Section  2.1.1.  The  rapidly  changing  vibration  dynamics  of  the  LV’s 
components  create  the  axial  and  lateral  vibration  loads.  The  axial  and  lateral  dynamics 
are  derived  in  Section  2.1.2.  Finally,  having  only  referred  to  the  environmental  loads  in 
general  terms,  they  are  fully  described  in  detail  in  Section  2.2. 

2.1.1  Gravity  Turn  Equations 

The  equations  of  the  trajectory  are  derived  in  the  form  of  a  gravity  turn  [26].  The 
gravity  turn  equations  are  needed  to  create  the  quasi-static  acceleration  loads  of  the 
launch  vehicle.  The  model  uses  a  non-rotating,  earth-centered  system.  For  computational 
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simplicity,  the  earth  is  assumed  to  be  spherical.  The  equations  for  this  study  are  derived 
for  two-dimensional  planar  motion.  The  assumed  symmetries  of  the  launch  vehicle  and 
satellite  allow  this  simplification. 

Figure  1  is  the  free-body  diagram  of  the  rocket  in  flight.  Ignoring  lift,  the  forces 
acting  on  the  LV  are  drag  (D),  thrust  (T)  and  weight  (Mg).  In  order  to  consider  the  LV  in 
the  familiar  local  horizon  frame  of  reference,  range  (X)  and  altitude  (H)  are  used  in  a 
non-rotating  earth  system.  Each  term  is  taken  along  one  axis  of  an  orthonormal  set  based 
at  the  earth’s  center.  The  range  direction,  X,  is  the  ei  component.  The  altitude,  H,  is  in  the 
e2  direction  and  the  e3  direction  is  the  axis  of  rotation  of  the  system  or  05.  Thus,  the 
instantaneous  position  of  the  LV  is  always  given  by  Equation  1.  Likewise,  the  angular 
velocity,  co,  is  given  in  Equation  2. 


v2 


Figure  1  Free  Body  Diagram  of  Launch  Vehicle  in  Flight 
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Taking  two  derivatives  of  the  position  vector  results  in  Equation  3.  This  result 
captures  the  centripetal  and  Coriolis  acceleration  of  the  rocket  as  it  flies  into  space.  When 
combined  with  the  thrust,  drag  and  weight  forces  acting  on  the  launch  vehicle  we  can 
analyze  the  trajectory  of  the  vehicle  and  also  compute  the  accelerations  on  the  payload. 
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The  flight  path  angle,  given  by  a,  indicates  the  direction  of  the  velocity  vector 
with  respect  to  the  local  horizon.  The  thrust  vector  and  drag  also  act  along  the  velocity 
vector.  Thus,  we  must  resolve  the  forces  acting  on  the  LV,  the  v  frame  of  reference,  along 
the  velocity  vector  to  simplify  the  equations  of  motion.  The  result  is  Equation  4. 
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Finally,  one  last  simplification  is  possible.  Noting  that  the  acceleration  can  be 
given  by  r  =  Vv,  +  Vocv2  in  the  LV  coordinate  frame,  the  equations  can  be  reduced  from 
two  second-order  equations  to  four  first-order  state  equations.  This  clearly  simplifies 
integration  of  this  highly  non-linear  system,  and  allows  access  to  meaningful  terms  like 
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the  vehicle’s  speed  and  axial  acceleration.  Equation  5  is  the  final  set  of  equations  in  body 
fixed  coordinates  (v  frame),  and  includes  the  necessary  integration  of  the  vehicle’s 
changing  mass. 
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The  equations  above  are  of  the  center  of  mass  of  the  LV.  The  forces  acting  on  the 
CM  are  characterized  below.  Realism  was  a  chief  goal  of  the  following  equations. 


Thrust  force  varies  for  two  reasons.  The  first  is  that  as  the  vehicle  ascends  into 
space  the  atmospheric  back  pressure  gradually  lessens.  This  has  the  effect  of  increasing 
thrust  force  at  higher  altitudes.  The  second  effect  is  that  solid  rocket  motors,  near  the  end 
of  stage  operation,  gradually  burn  out. 

Modeling  atmospheric  pressure  variability  is  accomplished  using  the  following 
relationship.  The  ambient  pressure  is  normalized  into  atmospheres  by  dividing  by  the 
constant  101,325  Pa. 


T  =Tsl  +(l-P/101325)(TVac  -Tsl)  (6) 

The  gradual  tailing-off  of  thrust  output  is  achieved  by  assuming  an  exponential 
decrease  in  thrust  during  the  last  seven  seconds  of  motor  operation.  This  models  a 
smoother  transition  to  a  zero  thrust  condition  at  the  end  of  the  bum.  Therefore,  the  source 
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of  shocks  at  staging  events  is  the  sudden  loss  of  mass  as  an  empty  stage  casing  is 
jettisoned  and  not  a  step  reduction  of  thrust.  The  tail-off  is  modeled  as  shown  in  Equation 
7.  Tto  is  the  reduced  thrust  (T  being  total  thrust)  which  begins  at  time  of  motor  cut-off, 
tco,  minus  seven  seconds. 

TT0=Te~W)]  (7) 


The  Drag  force  is  the  well  recognized  relationship  in  Equation  8  [21]. 


D  =  CD^pV2A 


(8) 


These  losses  in  performance  dramatically  reduce  the  payload  to  orbit  capability  of 
this  notional  launch  vehicle.  The  performance  of  the  notional  LV  modeled  here  is 
representative  of  actual  systems. 


2.1.2  Axial  and  Lateral  Vibration  Models 

The  next  step  is  to  create  a  vibration  model  representative  of  some  of  the  vibration 
modes  of  the  launch  vehicle  as  well  as  the  some  modes  of  the  satellite.  It  would  be  easy 
to  capture  as  many  LV  and  satellite  modes  as  possible,  but  computationally  costly.  Thus, 
very  early  on,  it  was  decided  to  limit  the  model  to  two  LV  modes,  one  satellite  mode  and 
the  isolator  critical  mode.  The  result  is  a  four  degree-of-freedom  (DoF)  system,  when  the 
rigid  body  mode  is  excluded.  The  lateral  system  will  be  similarly  separated. 

The  notional  launch  vehicle  model,  however,  is  of  a  staged  launch  vehicle  with 
the  potential  for  each  stage  to  have  its  own  characteristics.  For  modeling  simplicity  a 
standard  model  was  used  and  the  characteristics  of  each  stage  are  passed  into  the  model 
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as  the  vehicle  changes  configuration  throughout  its  flight.  The  model  that  results  is  a 
fixed  set  of  parameters  for  each  of  the  flight  configurations,  thus  allowing  continuous 
integration  as  the  vehicle  states  change.  This  has  the  nice  effect  of  minimizing  the 
number  of  simultaneously  running  integrators  and  reducing  the  number  of  blocks 
simulating  the  vehicle.  This  reduces  program  size  and  visual  complexity,  too. 

The  model  is  based  on  the  idea  that  the  operating  motor  or  engine  and 
aeroacoustic  noise  are  the  primary  sources  of  vibration  throughout  the  flight.  Therefore, 
the  LY  is  divided  into  three  parts  -  the  Lower  (L)  and  Upper  (U)  components  of  the 
operating  stage,  and  the  Inert  (I)  component  comprised  of  the  remaining  stages  of  the  LV 
that  have  not  been  burned.  These  divisions  allow  for  a  DoF  for  the  motor  vibration  source 
and  the  aeroacoustic  noise  source.  What  remains  of  the  model  is  the  satellite  and  the 
isolation  system.  The  satellite  is  modeled  as  a  mostly  solid  body  with  one  axially  moving 
part  and  one  laterally  rotating  part  of  some  specified  mass.  This  allows  the  modeler  to 
look  at  some  satellite  mode  of  interest.  The  isolation  system  (if  there  is  one)  is  modeled 
as  a  spring  or  a  spring  and  an  active  controller  system  separating  the  LV  and  the  satellite. 
The  satellite’s  degree  of  freedom  within  the  model  on  top  of  the  isolator  is  referred  to  as 
the  Payload  (P)  component  and  the  payload  sub-components  are  designated  as  Sub¬ 
component  1  (1)  in  the  axial  direction,  and  Sub-component  2  (2)  in  the  lateral  direction. 
Figure  2  demonstrates  the  transition  from  a  notional  launch  vehicle  design  to  the  standard 
model  concept.  The  payload  sub-components  are  not  visible  in  this  presentation  because 
they  are  within  the  satellite  component. 
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Figure  2  Notional  Launch  Vehicle  to  Vibration  Model 

2.1.2.1  The  Axial  Equations 

Based  on  the  conceptual  model  described  above  it  is  possible  to  derive  a  system 
of  equations  to  describe  the  motion  of  the  vehicle  components. 

The  axial  (vi)  equation  of  motion  (Equation  5)  describes  the  bulk  motion  of  the 
launch  vehicle.  Summing  forces  in  the  Vi  direction  on  each  of  the  LV  components  results 
in  the  five  equations  of  Equation  11.  For  brevity,  terms  such  as  KLu,  representing  the 
stiffness  force  between  the  L  and  U  components,  are  an  abbreviation  for  the  stiffness 
force  equation  kLU  (XL  -  Xv  ). 


2-7 


(9a) 


mi(X  l  ^  g  )  T  N  M  mL  < 


f  y  ^  ^ 

g - - — 

R*+H 


X  H  f 
Isin  + — - — cos  >-K,„  -C 


Rm+H 


LU  LU 


/  ^  2  A 


My(Xu  +  Xq  )  —  /Wy'll  g 


R~+H 


X  H  I 

lsin  d  °  zj  cos  I  —  ^ui  ~~  ^uj  +  Kuj  +  Cvu 
Am  H”  Jnl 


wij(X j  +  X  G )  —  ftij  < 


(  v  2  A 

g~ 


R*+H 


mp(X p  +  XG  )  =  -mF 


ml(Xl  +  XG)  =  -ml< 


'  XG>  ' 


xgh 

sin  H — — - cos 

Rn+H 

XrH 


X  ip  C  [/,  +  Kut  +  Cw  D  +  Na  I  l  Ir 


g~ 


r9+h 

2  A 


sin  +— — — — —  cos  \  —  K  F1  —  C  pi  +  K  ip  +  C  jp  + 1  ^  + 1 R 


Ra+H 


X 


Ra,+H 


X  H 

sin  H - - - cos  \  +  Km+CB 

Rm  +  H 


(9b) 

(9c) 

(9d) 

(9e) 


l/>  i  1  '-'pi 


Manipulating  the  equations  to  solve  for  the  component  accelerations,  and 
substituting  from  Equation  4  for  the  acceleration  of  the  vehicle  center  of  mass  yields  the 
following. 
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Dividing  each  equation  of  Equation  10  by  the  mass  of  each  component  solves  for 
the  accelerations.  This  system  of  equations  will  be  easily  programmed  into  Simulink. 
Next,  the  lateral  equations  of  motion  are  derived. 
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2.1.2.2  The  Lateral  Equations 


Following  a  similar  line  of  reasoning,  the  lateral  equations  of  motion  are  derived. 
The  equations  developed  here  will  allow  analysis  of  the  rocking  motion  of  the  satellite 
caused  by  vibration  and  transient  events  like  wind  gusts. 

Referring  again  to  Equation  2,  we  take  one  derivative  of  the  angular  velocity  and 
get  the  result  in  Equation  11.  Figure  3is  the  free  body  diagram  of  the  moments  acting  on 
the  launch  vehicle.  Multiplying  this  by  the  total  moment  of  inertia  of  the  notional  launch 
vehicle  and  adding  moments  caused  by  vibration  and  gust  forces  to  the  moment  equation 
yields  the  result  in  Equation  12. 


NmTn1n+G1n  v2 
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Figure  3  Free  Body  Diagram  of  Moments  Acting  on  the  LV 
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Analyzing  the  internal  moments  of  the  LY  is  similar  to  the  axial  force  analysis. 
Each  of  the  rocket  body  components  is  considered  with  respect  to  the  vehicle’s  center  of 
mass.  This  simplifies  the  final  equations  of  motions,  but  adds  the  requirement  for 
tracking  the  motion  of  the  center  of  mass  within  the  LV.  This  constraint  is  constantly 
changing  position  during  motor  operation.  Thus,  a  necessity  to  the  final  model  will  be  a 
full  description  of  the  rocket  components.  Equation  13  lists  the  lateral  equations. 

iA+cujeL-®,)+KU&L-^)=K  <13a> 

CMU1  ”  ®7  )  ^MUI  (®t/  —  ®  =  CMLU  (®L  +  ^ MLU  (®L  —  ^ 

/A  + 1  cIP(e,  -  eP)+^KIP(Q,  -  &P)  =  cMUI{@u  -  ©,) + kmui{©v  -  ©,) +m,  (13c) 

bQp  +  cMP2(eP  -  ©2) + kMP2(&p  -  ©2) = ^  c,P(©,  -  ePy +  ±  kip(&,  -  eP) + mp  (1 3d) 

i2e2  -  cP2(©p  -  02)  -  kp2(©p  -  02)  =  m2  (i3e) 

The  Mi  terms  are  the  internal  torques  and  are  expanded  in  equation  14. 
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Except  for  II  and  Ir  the  /,  terms  are  the  mass  moments  of  inertia  relative  to  the  LV 
center  of  mass  (c.m.).  The  /y  term  is  the  LV  total  mass  moment  of  inertia.  IL  and  IR  are 
the  left  and  right  lateral  control  force  inputs.  NM  is  the  vibration  noise  generated  by  the 
operating  motor,  and  Na  is  the  aeroacoustic  vibration.  The  motor  noise  is  observed  at  the 
nozzle  of  the  operating  stage  and  is  at  a  distance  In  from  the  c.m.  of  the  LV.  The  lateral 
vibration  is  assumed  to  be  a  fractional  value  of  the  axially  generated  vibration  and  is 
multiplied  by  the  lateral  transmission  factor,  Tn-  Similarly,  the  aeroacoustic  noise  is 
observed  at  the  payload  interface,  a  distance  //  from  the  LV  c.m.,  and  is  multiplied  by  the 
interface  lateral  transmission  factor,  Tj. 

In  several  places  in  Equations  13  and  14  terms  with  d/4  appear.  These  terms  relate 
to  the  interface,  where  d  is  the  diameter  of  the  PAF.  The  model  must  represent  a  circular 
distribution  of  the  interface  stiffness,  such  as  when  a  satellite  is  mounted  on  a  ring  of 
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passive  isolators  at  the  payload  interface.  The  derivation  for  this  is  rather  simple.  Figure 
4below  is  a  diagram  of  the  parts  of  the  interface. 


Figure  4  Diagram  for  the  derivation  of  the  lateral  stiffness 

The  derivation  is  as  follows.  Assuming  that  the  relative  deflection  angle,  a,  of  the 

payload  to  the  interface  is  small,  then  the  following  linearized  relationship  is  derived. 


dM  .  .  2 
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Equation  15  states  that  the  rate  of  change  of  the  moment  at  the  interface  with  respect  to 
the  incremental  angular  coordinate  on  the  payload  interface  ring  is  equal  to  the 
incremental  circumferential  stiffness,  k,  at  a  distance  rcos ©  from  the  center  of  rotation 
times  the  deflection  of  k  by  the  arc  of  arcos©.  Integrating  this  result  around  the  payload 
interface  ring  and  simplifying  yields  the  following  relationship  between  the  axial 
stiffness,  Kn>,  and  the  equivalent  stiffness  moment. 
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M  =  r 


(16) 


a 


It  is  also  necessary  to  develop  a  method  for  computing  the  moments  of  inertia  of 
the  components  of  the  LV.  For  the  solid  rocket  motor  the  most  common  design  is  a 
progressive  burn  design  similar  to  Figure  5.  This  calls  for  special  consideration  as  the 
motor  bums  and  the  mass  moment  of  inertia  is  changing  at  a  rate  proportional  to  the  mass 
flow  rate  of  the  operating  motor.  Also  there  are  mass  distribution  considerations  with  the 
inclusion  of  a  nozzle.  The  model  derived  also  combines  several  stages  into  one  degree  of 
freedom. 

Equation  17  was  used  to  develop  a  method  for  solving  the  instantaneous  mass 
moment  of  inertia  about  the  c.m.  of  the  component  [10].  The  equation  is  for  a  cylindrical 
solid.  Thus  a  shell  for  the  motor  case  has  to  be  calculated  and  then  a  solid  for  the 
propellant  added  back.  Likewise,  a  portion  of  the  propellant  has  to  be  removed 
incrementally  as  the  motor  burns.  Finally,  adding  the  mass  of  the  nozzle  contributes  to 
the  total  moment  of  inertia.  Equation  15  is  the  method  of  calculation. 
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I  Compnent  ^  CaseVolume 


I CaseVoid  ^Pr opeUuntVolume  ^  ?t  opellantVoid  *  Nozzle 


(18) 


For  each  component  the  volume  of  the  case,  the  volume  of  propellant  and  wall 
thickness,  and  nozzle  mass  must  be  specified  to  fully  define  the  component.  Also,  the 
center  of  mass  of  each  component  must  be  instantaneously  computed  to  properly 
calculate  the  contribution  of  the  nozzle. 


It  is  important  to  note  at  this  point  that  Equation  16  allows  the  decoupling  of  the 
axial  and  lateral  systems  of  equations.  They  are  merely  linked  by  the  axial  stiffness  of  the 
interface.  This  is  only  true  for  the  simplified  lumped  mass  model  developed  here.  If  the 
system  were  modeled  as  beams  in  bending  and  compression  then  there  would  be  a  direct 
coupling  of  lateral  stiffness  to  axial  compression  due  to  the  acceleration  of  the  LV  [7]. 
The  use  of  linear  springs  in  the  model  eliminates  this  effect. 


2.2  Environmental  Loads 

Three  types  of  loads,  other  than  the  quasi-static  loads  of  gravity,  thrust  and  drag 
are  considered.  They  are  motor  generated  vibration,  aeroacoustically  generated  vibration 
and  gust  loads.  Each  of  these  is  empirically  defined  so  that  a  representative  amount  of 
vibration  and  load  is  passed  to  the  satellite  at  the  PAF.  The  vibration  loads  presented  here 
are  filtered  for  frequency  content  in  Simulink.  A  fourth-order  Butterworth  high-pass  filter 
is  used  to  shape  the  frequency  output  of  a  band-limited  white  noise  source.  The  purpose 
is  to  more  closely  match  the  two-sigma  spectral  density  proposed  by  Larson  and  Newell 


for  a  range  of  launch  vehicles  [12].  Also,  the  frequencies  of  interest  are  limited  to  the  0- 
100  hertz  range  [11].  This  limitation  helps  to  reduce  the  stiffness  of  the  equations  to  be 
integrated. 

2.2.1  Motor  Generated  Vibration 

The  motor  generated  vibration  source  is  defined  to  originate  at  the  nozzle  of  the 
operating  stage.  It  includes  a  random,  band-limited  white  noise  component  and  a  resonant 
bum  component.  The  resonant  burn  component  is  assigned  a  frequency  and  amplitude 
that  is  representative  of  the  motor  being  considered.  The  band-limited  white  noise 
component  is  also  assigned  a  magnitude  that  results  in  an  appropriate  level  of  vibration  at 
the  payload  attach  fitting  (PAF).  The  random  component  is  also  weighted  by  a  ratio  of 
the  current  thrust  to  the  maximum  thrust  of  the  first  stage  motor.  Thus,  later  stages 
produce  less  noise  than  the  first  and  largest  stage.  This  is  true  in  actual  operation  [11]. 

Resonant  burn  is  a  condition  arising  in  solid  rocket  motors  when  the  combustion 
cavity  and  the  mass  flow  rate  create  a  resonance.  Resonant  bum  is  a  problem  for  large 
solid  rocket  motors.  The  vibrations  cause  an  oscillation  of  the  propellant  remaining  in  the 
motor  [11]. 

2.2.2  Aeroacoustically  Generated  Vibration 

For  the  purpose  of  this  model  the  atmosphere  is  defined  to  be  isothermal.  The 
temperature  is  based  on  the  1976  Standard  Atmosphere.  An  altitude  averaged  temperature 
was  computed  to  be  249K.  The  result  is  a  smooth  pressure  and  density  contour  into 
space.  Neither  value  will  actually  go  to  zero  but  this  is  true  for  the  LEO  orbits  considered 
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by  this  model.  Thus,  the  Mach  Number  of  the  LV  can  be  expressed  by  M=v/c.  Where  the 
speed  of  sound,  c,  is  constant  for  all  altitudes.  This  reduces  computational  demands  and 
fits  well  with  the  simplified  drag  model  also  in  use.  As  a  result  of  drag  the  model 
demonstrates  the  orbital  instability  caused  by  drag  when  allowed  to  run  for  long  times. 

Several  sources  site  the  transonic  flight  regime  as  being  the  most  intense  period  of 
aeroacoustic  loading  [9],  [11].  Another  period  of  intense  loading  is  maximum  dynamic 
pressure,  or  max-q.  To  that  end  the  vibration  load  acting  on  the  upper  most  part  of  the  LV 
structure  is  assigned  magnitude  according  to  a  profile  of  the  Mach  number  and  dynamic 
pressure.  This  relationship  is  given  as  inversely  proportional  to  Mach  number  and 
proportional  to  the  dynamic  pressure  and  density.  Empirical  tests  of  the  model  to  create 
appropriate  accelerations  at  the  payload  interface  dictated  the  constants  used.  Thus,  the 
aeroacoustic  vibration  scale  factor  is  given  as  Equation  19. 


S  =  p 


— +  |l-M2'05 


+  12- 


101325 


(19) 


This  equation  is  plotted  for  first  stage  flight  in  Figure  6.  Note  the  peak  at 
approximately  27  seconds.  It  corresponds  to  passage  through  Mach  1. 
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Aeroacoustic  Scale  Factor  vs  Time 


Figure  6  Aeroacoustic  Vibration  Scale  Factor  vs  Time 


2.2.3  Gust  Loads 

The  gust  load  for  the  lateral  case  is  modeled  as  a  step  input  with  duration  and 
magnitude  tuned  to  create  a  resonance  with  a  specified  frequency.  This  is  a  typical 
coupled  loads  analysis  approach  to  creating  a  gust  load  [1 1].  In  this  model  the  gust  load  is 
input  at  the  base  of  the  Lower  component. 

2.2.4  Model  Mission  Profile 

The  discussion  of  the  model  environment  would  be  incomplete  without 
characterizing  the  mission  profile.  The  profile  developed  for  this  study  is  similar  to  small 
solid-propellant  launch  vehicles.  Figure  7  shows  the  event  timing  of  the  LV  ascent.  The 
trajectory  modeled  is  an  eastern  equatorial  launch  so  that  planar  motion  can  be  assumed. 
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Figure  7  LV  mission  profile 

Further,  the  basic  vehicle  parameters  are  listed  in  Table  1.  These  characteristics 
were  required  to  accurately  model  the  LV. 

Combining  all  of  the  elements  discussed  in  this  chapter  results  in  the  total  mission 
acceleration  profile  that  appears  in  Figure  8.  This  acceleration  profile  was  generated  by 
the  model  developed  for  this  research. 
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First 

Second 

Third 

Fourth 

Payload 

PLF 

Length  (m) 

9 

8.9 

2.7 

1.23 

2 

6 

Diameter  (m) 

2.4 

1.3 

1.3 

1.3 

1.3 

1.3 

Total  Mass  (kg) 

53020 

14020 

3370 

985 

1780 

205 

Propellant  Mass  (kg) 

48809 

12152 

3025 

782 

Case  Mass  (kg) 

3415.7 

1631.9 

226.3 

150.2 

Nozzle  Mass  (kg) 

795.3 

235.4 

118.7 

52.8 

Mass  Flow  Rate  (kg  sec'1) 

-598.883 

-166.011 

-37.485 

-11.333 

Burn  Time  (sec) 

81.5 

73.2 

80.7 

69 

Initial  Grain  Dia.  (m) 

0.4 

0.3 

0.15 

0.1 

Propellant  Density  (kg  m‘3) 

1354.013 

1277.124 

960.723 

614.027 

Case  Density  (kg  m'3) 

8010.822 

7911.227 

8001.177 

7996.407 

Case  Wall  Thickness  (m) 

0.0063 

0.0057 

0.00257 

0.00375 

Table  1  LV  Physical  Properties 
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2.3  Controller  Theories 


The  controller  theories  studied  for  this  thesis  are  briefly  discussed.  Proportional 
(P),  proportional  plus  derivative  (PD),  proportional  plus  integral  (PI),  and  proportional 
plus  integral  plus  derivative  (PID)  are  the  classic  control  methods  explored  as  alternatives 
for  hybrid  isolation.  Positive  position  feedback  (PPF  was  also  studied. 

2.3.1  PI,  PD,  and  PID 

PI,  PD  ad  PID  are  the  standards  of  discussion  for  most  classical  control 
approaches.  These  controllers  are  developed  for  negative  closed-loop  feedback.  Many  of 
the  controller  design  methods  are  prescribed  for  this  controller  family.  Tools  like  root- 
locus  and  Bode  diagrams  have  detailed  design  methodologies  to  support  PID  control. 

The  fundamental  PID  controller  feedback  equation  is  shown  below  in  Equation 
20.  Proportional  control  retains  only  the  gain  K.  The  Kis'1  is  the  integral  term,  and  the 
term  Kds  is  the  derivative  control  term. 

G(s)  =  K  +  ^-  +  Kds  (20) 

■S' 

For  a  more  detailed  discussion  of  PI,  PD,  and  PID  control  and  other  classical 
control  methods,  please  refer  to  Reference  [5]. 

2.3.2  Positive  Position  Feedback 

The  positive  position  feedback  (PPF)  control  law  is  shown  in  Equation  21.  The 
first  equation  listed  is  the  equation  of  motion  of  the  system.  The  disturbance  force  is  d, 
and  control  force  is  yo)n2Z.  The  active  feedback  equation  specifies  the  frequency  to  be 
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isolated,  coc,  and  takes  its  input  from  the  position,  X,  of  the  payload  component.  The 


appearance  of  the  PPF  equations  is  similar  to  a  2DoF  system,  however  Equation  21b  is  of 
a  virtual  mass. 

X  +  2C«„  X+a>]X=  ya>lZ + d  (21a) 

Z  +  2£c0)cZ  +  CQ2nZ=CQ*X  (21b) 

For  a  more  detailed  discussion  of  PPF  control,  please  refer  to  Goh  and  Caughey 
(1985)  [6], 

2.4  Model  Layout 

The  top-level  Simulink  model  is  briefly  addressed  in  this  section.  The  model  is 
discussed  in  nauseating  detail  in  Appendix  A.  Appendix  B  provides  the  source  code 
listing  of  the  various  models  implemented  for  the  analysis  presented  in  Section  3. 


model.  The  clock  drives  many  events  of  the  model.  Staging,  gusts  and  other  flight  events 
are  triggered  by  the  clock.  The  model  divides  roughly  into  two  parts  -  environment 
generation  and  dynamic  output.  The  Noise  Generator,  LV  CG  Dynamics,  and  Mass 
Properties  blocks  generate  the  environment  that  excites  the  model  vibration  dynamics. 
The  Axial  Forcing  Functions  and  the  Axial  Vibration  Dynamics  calculate  the  dynamic 
output  of  the  model.  The  end  results  are  the  positions,  velocities  and  accelerations  of  the 
LV  components.  The  model  demonstrates  this  flow  from  clock  to  dynamic  output.  The 
blocks  of  the  model  can  be  combined  to  create  new  models. 

Environment  generation  in  any  simulation  is  the  same.  The  three  blocks,  LV  CG 
Dynamics,  Mass  Properties  and  the  Noise  Generator,  are  used  universally.  These  blocks 
represent  the  characteristics  of  the  launch  vehicle  and  its  mission  profile.  These  are  the 
only  three  time  dependent  blocks  and  are  the  source  of  the  full  simulation’s  non-linearity. 
These  will  be  discussed  later  in  Appendix  A.  The  LV  CG  Dynamics  block  is  based  on  the 
gravity-turn  equations  derived  in  Section  2.1.1. 

The  most  basic  model  is  the  five-DoF  simulation  in  Figure  9.  Adding  active 
controls  results  in  models  similar  to  Figure  10.  This  model  uses  model  reference  adaptive 
control.  The  LV  model  can  be  radically  simplified  to  use  a  state-space  linear  time 
invariant  representation  for  rapid  controller  prototyping  as  in  Figure  1 1 .  Note  that  in  the 
state-space  model  the  dynamic  output  does  not  rely  on  the  different  component  masses 
due  to  time-invariance.  This  model  approximates  the  vehicle’s  state  at  27  seconds  after 
launch.  The  state-space  block  calls  on  workspace  variables  defined  in  sscont3ppf.m 
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(Appendix  B)  executed  prior  to  the  simulation  that  establish  the  time  invariant  model 
components. 


Axial  Forcing 


Figure  11  Estimator  model  of  hybrid  isolation  system  using  PPF  control 


2.5  Summary  of  Axial  and  Lateral  Equations 

The  equations  of  lateral  and  axial  motion  compared  side-by-side.  The  reader  will 
note  that  the  only  dependence  the  lateral  equations  have  on  the  axial  equations  is  the 
stiffness  of  the  payload  interface.  This  fact  is  substantiated  by  the  assumption  of  a  linear 
spring  separating  the  Payload  and  Inert  model  components. 
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(10b) 

(10c) 

(lOd) 

(lOe) 


A  =  KP ,  +Cpl - l-(T  +  Nm  -D  +  Na) 

mG 

Equation  10  above  is  the  axial  equations  and  Equation  13  below  is  the  lateral 
equations  of  motion. 


/A  =  ~CMLU  (®l-®«)-kmlu(ql-qu)+ml 

^V®V  ~  ~^MV1  (®t 1  ~  ®/ )  —  ^  MU!  (®u  ~  ®/  )  CMLu(@L  —  Gy  )  +  Kmlv  (0l  —  0y  )  +  My 
I  i®i=~Cip(@1-Qp)-^Kip(®i-®p)  +  Cmui(Qu-®i)  +  Kmui(Qu-@i)  +  Mi 

ipgp=-cMP2(Qp -g2)-kMP2(gp -e2)+^cIP(@, -gp)+^kip(g, -gp)+mp 
I2Q2=Cp2(Qp -®2)  +  Kp2(®p  -02)  +  M2 


(13a) 

(13b) 

(13c) 

(13d) 

(13e) 
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3  Isolation  Results 


The  results  and  analysis  of  data  collected  using  the  models  described  in  the 
previous  section  is  presented.  The  aim  of  the  research  was  to  show  that  passive  and 
hybrid  isolation  of  whole  spacecraft  from  the  LV  induced  vibration  environment 
substantially  reduce  the  loads  felt  by  the  spacecraft  and  a  payload  sub-component.  While 
equations  of  the  axial  and  lateral  case  were  derived,  it  was  shown  that  the  axial  and 
lateral  cases  de-couple.  Thus,  the  focus  of  this  paper  will  be  on  the  axial  case  alone. 

This  section  is  organized  in  the  following  fashion.  After  developing  the  baseline 
response  of  the  model  (Section  3.1),  the  passive  isolation  results  will  be  analyzed 
(Section  3.2).  The  hybrid  results  will  then  be  presented  in  Section  3.3.  The  isolation 
system  effectiveness  is  compared  to  the  baseline  system  repeatedly.  In  each  case  the 
analysis  will  focus  on  the  first  90  seconds  of  flight.  This  time  period  captures  ignition, 
transonic  flight,  max-q,  and  staging/ignition  of  stage  two.  In  addition,  the  chief  variable 
in  the  performance  analysis  is  the  passive  isolator  critical  frequency.  The  analysis  of  the 
model  will  focus  on  passive  isolator  frequencies  between  25  and  40  Hz.  This  band 
brackets  the  second  mode  range  during  the  transonic  and  max-q  flight  periods. 

Several  methods  of  comparison  will  be  used  to  analyze  the  effectiveness  of  a 
particular  isolation  combination.  Power  spectral  density  (PSD)  comparisons  of  the 
acceleration  of  the  Baseline,  Passive  and  Hybrid  systems  for  the  Payload  and  Sub¬ 
component  are  used.  Supplementing  the  PSD  is  the  spectrogram.  The  spectrogram  is  a 
time  history  of  the  PSD  of  the  model.  It  takes  small  time  slices  of  the  output  signal  and 
performs  a  fast  Fourier  transform  of  those  slices  to  create  a  plot  mapping  the  intensity  of 
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a  given  frequency  band  with  respect  to  time.  Finally,  the  root-mean-sum-square  (RMS)  of 
the  acceleration  of  a  given  component  is  computed  for  each  run.  Comparison  of  the  RMS 
boils  the  performance  of  an  isolator  combination  down  to  a  single  number. 

3.1  Baseline  Characterization 

Of  the  five  degrees  of  freedom  of  the  model,  the  spacecraft  comprises  two.  The 
satellite  modeled  has  a  mass  of  1780  kg.  The  Sub-component  of  the  payload  is  1  kg.  The 
Sub-component  has  as  its  critical  vibration  frequency,  50  Hz.  This  was  chosen  to  match 
the  resonant  bum  frequency  (also  50  Hz)  of  the  first  stage. 

Frequency  vs  Time 


Given  the  non-linearity  and  time-variance  of  the  model,  some  expectation  of  the 
behavior  of  the  model  is  needed.  The  mass  of  the  LV  changes  constantly  as  the  motor 
operates.  This  causes  the  mass  of  the  Lower  and  Upper  component  of  the  model  to 
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lessen.  This  lessening  of  mass  combined  with  constant  structural  stiffness  of  the  LV  will 
cause  the  LV  modes  to  constantly  vary  throughout  stage  operation.  It  will  be  necessary  to 
consider  this  variation  and  the  multiple  modes  of  the  LV  structure  to  select  a  successful 
isolation  design.  Figure  12  plots  the  time-varying  modes  of  the  LV.  It  supports  the 
previous  conclusion  and  the  results  of  past  analysis  of  launch  vehicle  structural  modes 
[7]. 

The  inter-component  stiffness  were  kuj  =  4xl08  N/m  and  kui  =  1x10s  N/m.  As  the 
diagram  above  shows  this  creates  a  low  frequency  first  mode  that  varies  only  slightly, 
and  a  rapidly  changing  higher  frequency  mode.  This  second  mode  has  the  important 
feature  of  transitioning  through  the  critical  frequency  of  the  model  Sub-component.  This 
combined  with  resonant  burn  of  the  first  stage  motor  creates  several  resonance  situations. 

The  baseline  LV/Payload  system  is  taken  from  the  five  DoF  model.  The  baseline 
is  intended  to  represent  a  spacecraft  that  is  hard-mounted  to  the  payload  attach  fitting 
(PAF).  The  PAF  is  generally  a  ring  to  which  the  spacecraft  separation  system  is  attached. 
The  payload  isolation  system  mounts  between  the  PAF  and  the  spacecraft  separation 
system. 

In  the  model,  the  hard-mounted  case  is  considered  to  be  an  extremely  stiff  linear 
spring  between  the  Inert  and  Payload  components.  The  critical  frequency  of  the  hard¬ 
mount  is  200  Hz.  The  first  figure  presented  is  the  spectrogram  of  the  baseline  system  (See 
Figure  13).  Clearly  visible  in  the  spectrogram  are  the  first  and  second  modes  of  the  LV, 
and  the  resonant  bum  condition,  which  exists  throughout  first  stage  operation.  At  80-82 
seconds  after  launch  the  second  stage  event  is  obvious  as  a  broadband  shock. 
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Intensification  at  approximately  60  seconds  after  launch  is  visible  along  the  50  Hz  line  as 
the  resonant  burn  and  second  LV  mode  cross.  The  black  band  along  the  right  edge  of  the 
graph  is  an  artifact  of  scaling  and  FFT  windowing. 

Spectrogram  of  Upper  Component 
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Figure  13  Baseline  Spectrogram  of  Upper  Component 

The  next  three  diagrams  are  power  spectral  density  (PSD)  plots  of  the  frequency 

response  of  the  Inert  component.  Payload  component  and  Sub-component.  The  PSDs  are 
time  averaged  frequency  response  plots.  It  is  easy  to  locate  the  LV  modes  and  in  the  case 
of  (Fig  SC)  the  critical  frequency  of  the  Sub-component  is  the  prominent  feature  of  the 
diagram. 
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Baseline  PSD  of  Inert  Component 
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Figure  14  Baseline  PSD  of  Inert  Component 

Baseline  PSD  of  Payload  Component 
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The  RMS  acceleration  in  g  of  the  Payload  and  Sub-component  of  the  baseline 
case  are  listed  in  Table  2. 


Baseline-PL 

Baseline-SC 

Acceleration  (g)|  13.35 

0.7716 

Table  2  Baseline  RMS  -  g  acceleration  of  the  Payload  and  Sub-component 

Spectrograms  of  the  Payload  and  Sub-component  indicate  the  environmental 

factors  that  most  influence  the  system  dynamics.  Figure  17  is  the  spectrogram  of  the 
baseline  system  for  the  Payload  component.  The  most  prominent  features  of  this  figure 
are  the  trace  of  the  first  LV  mode,  the  intensification  of  the  50  Hz  excitation  at 
approximately  50  seconds,  and  the  broadband  impact  of  the  staging  event  at  81.5 
seconds. 
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Spectrogram  of  Payload  Component 
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Figure  17  Baseline  Spectrogram  of  Payload  Component 


The  baseline  spectrogram  of  the  Sub-component  is  included  in  Figure  18. 


Spectrogram  of  Payload  Sub- component 
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Figure  18  Baseline  Spectrogram  of  the  Sub-component 
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The  most  significant  feature  of  the  spectrogram  of  the  Sub-component  is  the 


appearance  of  a  mode  at  7  Hz.  This  mode  was  not  predicted  by  analysis  of  the  system. 


3.2  Passive  Isolation  Results 

Use  of  passive  isolation  in  the  full  model  is  marked  by  significant  load  reduction 
for  the  whole  spacecraft.  The  Sub-component  experiences  a  minor  increase  in  loading, 
but  for  the  most  part  remains  near  the  un-isolated  response  level. 


3.2.1  Overall  System  Response 

The  passive  isolator  was  studied  for  first  stage  flight,  from  zero  to  90  seconds.  A 
range  of  passive  isolator  frequencies  (Critical  Frequencies:  25-40  Hz)  was  studied  so  that 
a  broader  understanding  of  the  effects  of  passive  isolation  was  possible.  Figure  19  shows 
the  RMS  -  g  acceleration  of  the  payload  and  Sub-component.  The  Payload  (or  whole 
spacecraft)  greatly  benefits  from  the  passive  isolator  at  all  frequencies  observed. 
However,  the  Sub-component  experiences  increasing  loading  as  the  isolator  frequency  is 
increased.  This  is  expected,  because  harmonic  excitation  from  the  resonant  burn  effect  is 
enhanced  by  the  resonance  of  the  passive  isolator. 
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RMS  Acceleration  vs  Passive  Isolator  Frequency 


—♦—Passive  Payload 
■hi —  Passive  Sub-component 
••■■A— Baseline  Payload 
•  x  Baseline  Sub-Component 


Figure  19  Passive  isolator  effectiveness  vs  passive  isolator  frequency 

Data  was  not  taken  in  the  vicinity  of  the  first  LV  mode  (13.2  Hz  at  lift-off), 

because  the  lateral  mode  of  the  Payload  component  with  a  low  frequency  passive  isolator 
might  be  low  enough  to  interact  with  the  guidance,  navigation  and  control  of  a  real  LV. 
Thus,  25  Hz  was  chosen  as  a  reasonable  lower  limit.  The  data  show  that  the  passive 
isolator  significantly  reduces  payload  vibration  for  the  entire  range  of  passive  isolator 
critical  frequencies.  In  addition,  there  is  a  marginal  increase  for  the  Sub-component. 
Figure  20  compares  the  passive  isolated  Sub-component  to  the  baseline  case.  The 
gradually  reduced  loading  continues  beyond  40  Hz.  At  70  Hz  the  Sub-component 
acceleration  is  less  than  that  of  the  baseline  case.  Figure  21  explains  this  result.  The 
crossover  of  the  Payload  baseline  and  passive  isolated  case  occurs  at  approximately 
60Hz.  Thus,  in  the  case  of  this  model,  satellite  components  with  resonance  above  60  Hz 
will  be  further  isolated,  but  below  60  Hz  satellite  sub-components  will  experience 
increased  loading. 
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RMS  Acceleration  of  the  Sub-component  vs  Passive  Isolator  Frequency 


Figure  20  Passive  Isolator  Effect  on  the  Sub-Component 


PSD  of  Passively  Isolated  Payload  Component 


Figure  21  PSD  indicating  60  Hz  crossover  frequency  on  the  Payload  Component 
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3.2.2  Best  Case  Response:  30  Hz  Passive  Isolator  Critical  Frequency 

The  best  result  is  obtained  at  a  passive  isolator  critical  frequency  of  30  Hz.  At  this 
frequency  the  passive  isolator  reduces  the  vibration  load  on  the  payload  87.5  percent 
from  the  baseline  case.  This  isolation  frequency  is  discussed  in  more  detail. 

Figure  22  is  the  PSD  of  the  Payload  component.  The  plot  compares  the  baseline 
and  passive  isolation  case  at  an  isolator  resonant  frequency  of  30  Hz.  The  passive  isolator 
effectively  attenuates  all  frequencies  above  30  Hz.  The  isolator  notably  reduces  the 
response  to  the  resonant  bum  excitation  at  50  Hz.  Figure  23  of  the  Sub-component 
further  confirms  the  earlier  discussion  that  the  passive  isolation  does  little  to  change  the 
response  of  the  Sub-component. 

PSD  of  Passively  Isolated  Payload  Component 
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The  spectrogram  of  the  Payload  component  in  Figure  24  clearly  shows  the  passive 
isolator  critical  frequency  at  30  Hz.  Comparisons  of  this  to  the  baseline  spectrogram  may 
seem  to  indicate  overall  intensified  loading  the  isolator  significantly  attenuates  high 
frequency  vibration  loading.  This  reduces  loading  by  87.5%  from  the  baseline  case  for 
the  Payload  component.  Figure  25  is  the  spectrogram  of  the  Sub-component.  The  passive 
isolator’s  effect  is  apparent  by  the  overall  leveling  of  the  loads  on  the  Sub-component. 

Spectrogram  of  Payload  Sub-component 
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Figure  25  Spectrogram  of  the  Passive  isolated  Sub-component 

3.2.3  Damping  Ratio  Effects  on  the  Passive  Isolator 

The  effect  of  damping  ratio  was  also  explored.  The  damping  ratio  was  allowed  to 
vary  between  0.25%  and  4%.  The  range  reflects  the  use  of  visco-elastic  damping 
applications  to  the  passive  isolation  system.  At  each  passive  isolator  critical  frequency 
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the  damping  ratio  was  doubled.  0  shows  that  over  the  available  range  the  damping  ratio 
made  little  difference  to  the  response  of  the  payload. 


3.3  Hybrid  Isolation  Results 

Extending  on  the  passive  isolation  concept  is  hybrid  isolation.  Hybrid  isolation  is 
the  combination  of  passive  and  active  vibration  control.  Most  hybrid  concepts  collocate 
the  passive  and  active  controllers.  One  suggested  design  is  the  placement  of  piezo¬ 
ceramic  actuator  stacks  in  the  high  strain  areas  of  the  passive  isolator  [19].  In  this  way, 
the  hybrid  design  seeks  to  maintain  the  fail-safe  capability  of  passive  isolation,  should  the 
active  component  fail. 
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Implementing  an  active  control  design  involves  taking  acceleration  data  with 
accelerometers  or  position  data  from  strain  calculations  based  on  output  from  the  piezo¬ 
ceramic  stacks.  This  data  is  in  turn  used  as  input  to  the  controller  scheme.  The  output  of 
the  controller  is  then  used  to  drive  the  active  component  of  the  isolation  system 

Again,  data  on  various  hybrid  designs  was  taken  across  a  range  of  critical 
frequencies  for  the  passive  isolator  (25-40  Hz).  This  variability  was  needed  to  explore  the 
synergy  of  passive  and  active  control.  At  each  passive  isolator  frequency,  the  active 
component  was  tuned  to  maximize  effectiveness.  Among  time-invariant  control  designs, 
proportional-plus-integral-plus-derivative  (PID)  and  positive  position  feedback  (PPF) 
were  attempted. 


3.3.1  Overall  System  Response 

An  interesting  result  of  the  hybrid  control  is  an  apparent  dependence  on  the 
passive  isolator  critical  frequency  for  the  amount  of  load  attenuation  that  is  gained  by 
hybrid  isolation.  At  some  frequencies  there  is  increased  loading  compared  to  the  passive 
isolation  case.  The  same  effect  was  observed  for  the  passive  case.  Figure  27  compares  the 
passive  isolated  and  hybrid  isolated  Payload  component.  Notably  the  passive  isolator 
outperforms  the  hybrid  isolator  and  conversely  the  hybrid  creates  gains  over  the  passive 
case 
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RMS  Acceleration  (g) 


Passive  vs  Hybrid  Isolation 


Passive-Payload 
~*~PPF-Payload 
. ^•••"PID-Payload 


Figure  27  .Passive  vs  Hybrid  Isolation  at  the  Payload  Component 


Over  the  range  of  passive  isolator  critical  frequencies  studied,  the  average 
improvement  was  an  83.7%  reduction  in  g-loading  on  the  whole  satellite.  Hybrid 
isolation  improves  on  passive  another  9.8%.As  was  mentioned  previously,  PID  and  PPF 
were  studied  as  potential  active  control  laws.  Figure  28  demonstrates  there  is  little 
difference  in  the  effectiveness  of  either  control  design.  Other  factors  outside  of  the  scope 
of  this  thesis  that  may  realistically  influence  control  law  choice  are  control  power 
required  to  effectively  deploy  the  design,  and  ease  of  implementation.  By  ease  of 
implementation,  recognize  that  PID  requires  accelerometers  installed  on  the  satellite, 
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while  PPF  will  require  the  use  of  strain  measurements  to  estimate  the  relative  position  of 
the  Payload  to  the  interface. 


PPF  vs  PID  Hybrid  Isolation 
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Figure  28  PPF  vs  PID  Hybrid  Isolation 

Figure  28  also  indicates  a  potential  risk  of  using  hybrid  isolation.  While 


tremendous  improvements  are  recognized  by  the  use  of  passive  isolation,  and  on  average 


hybrid  isolation  increases  attenuation  on  the  Payload,  the  Sub-component  is  exposed  to  a 
greater  level  of  loading.  This  effect  is  difficult  to  explain.  A  potential  explanation  may  be 


that  the  active  component  creates  harmonic  excitation  at  a  frequency  that  is  a  multiple  of 


the  critical  frequency  of  the  Sub-component.  Another  cause  may  be  that  the  active 


controller  creates  a  steady-state  error  that  is  transmitted  to  the  Sub-component. 
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3.3.2  Best  Case  Response:  32  Hz  Passive  Isolator  Critical  Frequency 

For  PPF  control,  the  next  two  plots  bear  this  last  concept  out.  Figure  29  compares 
the  ratio  of  the  PSD  of  the  passively  isolated  and  the  hybrid  (PPF)  isolated  Payload 
component  to  the  baseline  un-isolated  Payload  component.  The  diagram  shows  the 
steady-state  error  introduced  by  the  hybrid  isolator,  and  it  also  shows  where  the  PPF 
controller  has  eliminated  the  response  of  the  first  LV  mode  and  has  reduced  the  response 
at  the  passive  isolator  critical  frequency.  The  passive  isolator  critical  frequency  in  this 
plot  is  32  Hz. 

PSD  of  Transmission  Ratio  of  Payload  Component 


Figure  29  Transmission  ratio  comparison:  Passive  and  Hybrid  (PPF)  Systems 

A  further  comparison  of  the  hybrid  and  passive  case  is  possible  by  examining  the 

ratio  of  the  PSDs  of  the  hybrid  and  passive.  This  is  plotted  in  Figure.  The  PPF  controller 
is  tuned  at  three  frequencies  -  the  first  LV  mode  (~13.2  Hz),  the  passive  isolator  critical 
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frequency  (-31.5  Hz)  and  the  resonant  frequency  of  the  Sub-component  (-49.5  Hz).  The 
plot  shows  that  the  active  controller  affects  the  frequency  response  in  the  vicinity  of  the 
frequencies  addressed  by  the  PPF  control  law. 


PSD  of  Transmission  Ratio  of  Payload  Component 


While  the  PPF  controller  introduces  substantial  steady-state  error,  the  following 
plots  of  the  PSDs  of  a  hybrid  system  using  PID  control  indicate  that  PID  control  adds  no 
steady-state  error  into  the  system.  Continuing  the  comparison  at  a  passive  isolator  critical 
frequency  of  32  Hz,  Figure  31  plots  the  transmission  ratios  of  the  passive  and  hybrid 
systems.  Unlike  Figure  29,  there  is  no  steady-state  error  in  the  system,  yet  PID  eliminates 
the  first  LV  mode,  like  the  PPF  controller.  Thus,  steady-state  error  introduced  by  the 
controller  cannot  be  the  culprit  for  the  increased  loading  on  the  Sub-component. 
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PSD  of  Transmission  Ratio  of  Payload  Component 


Frequency  (Hz) 


Figure  32  clearly  shows  how  the  PID  controller  notches  the  first  LV  mode  (13.2 
Hz)  and  by  lowering  the  effective  frequency  of  the  passive  isolator  creates  a  symmetrical 
amplification  then  attenuation  around  the  passive  isolator.  The  gains  to  be  gotten  from 
PID  will  be  from  how  aggressively  the  controller  attenuates  the  first  LV  mode. 


Ratio  of  Spectrograms  of  PID  to  PPF 
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Figure  33  Ratio  of  the  Spectrograms  of  PID  to  PPF  at  32  Hz 

The  ratio  of  the  spectrograms  of  the  PID  controlled  isolator  to  the  PPF  controlled 

isolator  indicates  the  differences  of  the  two  control  methods.  Interpreting  the  spectrogram 
ratio  is  easy.  Areas  of  brightening  on  the  plot  indicate  that  the  PPF  controller  is 
outperforming  PID.  Conversely,  darkening  indicates  the  PID  method  is  outperforming 
PPF.  The  steady-state  error  induced  by  PPF  is  highlighted  along  the  bottom  edge  of  the 
plot. 
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3.4  Isolation  Summary 


Passive  isolation  is  characterized  by  a  significant  reduction  in  loading  on  the 
Payload.  The  Sub-component  experiences  a  marginal  increase  in  acceleration  loads,  but 
these  loads  decrease  with  the  critical  frequency  of  the  Sub-component.  Above 
approximately  60  Hz,  the  Sub-component  RMS  acceleration  drops  below  the  baseline 
acceleration  of  the  Sub-component.  This  might  encourage  designers  to  produce 
components  that  have  initial  modes  above  the  crossover  frequency  of  the  satellite  system 
in  question.  This  is  impractical  especially  in  light  of  one  of  the  fundamental  goals  of 
isolation  research  —  the  creation  a  launch  environment  that  encourages  standard  off  the 
shelf  components  for  use  on  any  satellite  system. 

Hybrid  isolation  promises  great  things  for  whole  satellite  isolation  by  further 
improvement  beyond  passive  isolation.  This  is  tempered  by  the  increased  response  of  the 
Sub-component  above  the  baseline  and  the  passive  cases.  Further  exploration  of  the 
modal  response  of  the  Sub-component  is  necessary  to  understand  this  phenomena. 
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4  Conclusions  and  Recommendations 


4.1  Conclusion 

This  study  follows  in  the  steps  of  the  recent  successful  demonstration  of  passive 
vibration  isolation  during  launch  of  a  satellite  payload.  Isolating  components  of  a  satellite 
payload  and  reinforcing  the  design  of  a  satellite  and  its  sub-components  has  long  been 
common  practice  in  the  satellite  manufacturing  community.  As  a  result,  much  study  has 
gone  into  improving  satellite  component  testing  and  satellite  construction  techniques. 
However,  as  satellite  payloads  become  smaller  in  order  to  take  advantage  of  cheap,  light 
capability  launch  vehicles,  the  susceptibility  of  these  small  satellites  to  launch  vibration 
loads  increases.  This  has  sparked  a  research  shift  in  the  vibration  isolation  community  to 
whole  spacecraft  isolation. 

The  research  of  the  past  several  years  in  launch  vibration  isolation  falls  roughly 
into  two  classes  of  modeling.  They  are  Coupled  Loads  Analysis  and  lumped  mass 
systems.  The  CLA  approach  sacrifices  speed  and  simplicity  for  accuracy.  While  the 
lumped  mass  approach  sacrifices  the  full  fidelity  of  CLA  for  broader  applicability  to 
general  design  considerations.  For  this  research  the  author  has  constructed  a  lumped  mass 
parameter  model  of  a  launch  vehicle  and  satellite  using  Simulink.  The  use  of  which  is  to 
construct  various  passive  and  hybrid  vibration  isolation  systems  and  study  their 
effectiveness. 

The  chief  difference  in  the  researcher’s  simulations  and  other  parallel  efforts  is 
that  the  model  constructed  is  of  a  full  LV  mission  profile.  CLA  and  other  lumped  mass 
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parameter  research  efforts  focus  on  very  short  simulations  of  less  than  20  seconds  of 
simulated  flight  around  very  specific  flight  environmental  conditions  or  in  many  cases 
much  simpler  loads,  such  as  impulse,  step  or  sine  sweeps.  Using  the  researcher’s  full 
flight  model  that  simulates  vibration  loads  encountered  during  launch  provides  insight  to 
the  demanding  environmental  transitions  for  which  a  hybrid  control  design  will  have  to 
account. 

Further,  the  model  constructed  for  this  study  also  captures  the  non-linear,  time- 
variant  nature  of  the  dynamics  of  a  LV  flight  to  orbit.  A  full  set  of  gravity  turn  equations 
is  the  basis  for  simulating  the  motion  of  a  LV  to  orbit  around  a  spherical  non-rotating 
earth.  The  changing  mass  properties  of  a  LV  are  represented  giving  rise  to  structural 
frequency  shifts  that  greatly  complicate  passive  and  hybrid  vibration  isolation  design. 
The  researcher  has  demonstrated  the  controllability  of  the  vibration  loads  on  a  satellite 
and  its  payload  component  for  axial  load  cases. 

Model  analysis  began  by  characterizing  a  baseline  case  against  which  the 
effectiveness  of  passive  and  hybrid  vibration  isolator  designs  could  be  compared. 
Subsequent  runs  of  passive  isolators  showed  tremendous  gains  over  the  un-isolated 
baseline  case.  The  improvement  was  on  average  an  80  percent  or  more  reduction  on  the 
whole  satellite.  Other  researchers  have  confirmed  this  level  of  improvement  [13],  [20], 
[27],  [28].  On  the  other  hand,  this  simulation  also  included  the  dynamics  of  a  sub¬ 
component  within  the  satellite.  The  surprising  result  was  the  dynamic  amplification  of  the 
loading  on  the  Sub-component.  Hybrid  isolation  designs  analyzed  by  the  researcher  PID 
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and  PPF  control.  The  research  indicate  that  PID  may  be  an  appropriate  control  technique 
for  whole  spacecraft  isolation. 

4.2  Recommendations  for  future  research 

The  relative  simplicity  of  the  five  DoF  model  designed  for  this  research  was 
sufficient  for  these  initial  results,  however  the  researcher  feels  that  future  work  with  this 
simulation  system  should  begin  by  increasing  the  model’s  degrees  of  freedom.  Doing  so 
will  allow  the  identification  of  more  LV  modes  and  satellite  modes.  The  first  stage  Inert 
component  of  this  model  is  very  large  and  should  be  broken  up  in  to  smaller  degrees  of 
freedom.  In  addition,  the  inclusion  of  additional,  and/or  off-centerline  axial  components 
in  a  lateral  model  will  create  interesting  coupling  effects  that  should  be  studied. 

While  this  study  explored  the  effect  of  varying  the  passive  isolator  critical 
frequency  it  does  not  analyze  the  conditions  that  arise  as  the  Payload  mass  or  Sub¬ 
component  mass  and/or  stiffness  change.  Collecting  this  information  will  characterize  the 
envelope  of  controllability  for  a  particular  design.  This  has  direct  ramifications  on  actual 
isolation  system  design  for  use  on  LVs.  Furthermore,  it  has  been  shown  that  the  interface 
stiffness  of  the  LV  has  a  direct  impact  on  the  control  force  that  can  be  generated  to 
reduce  vibration  loads. 

Adaptive  control  methods  should  be  explored.  Several  other  self-tuning  control 
methods  exist.  The  simulation  model  created  for  this  study  is  fully  modular  and  can 
easily  accept  new  control  designs.  The  only  limitation  to  what  can  be  tried  is  the  future 
researcher’s  creativity. 
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The  lateral  equations  of  motion  were  presented  in  this  study  as  part  of  the 
derivation.  This  model  focused  on  axial  vibration  isolation  once  it  was  shown  that,  for  the 
assumption  of  linear  spring  rates  between  the  LV  and  satellite  components,  the  axial  and 
lateral  vibration  cases  de-couple.  The  GFO  launch,  while  successfully  demonstrating 
passive  axial  vibration  isolation,  showed  that  lateral  loads  increase  as  axial  stiffness 
decreases  due  to  the  passive  isolator.  Laterally  stiffening  the  satellite  on  the  isolation 
system,  but  axially  softening  the  ride  will  limit  this  increased  rocking  mode.  Too  much 
lateral  motion  can  result  in  the  catastrophic  loss  of  the  satellite  if  it  were  to  strike  the 
payload  fairing. 

Identifying  the  source  of  the  dynamic  amplification  of  loads  on  the  Sub¬ 
component  caused  by  the  hybrid  isolator  is  a  priority.  Once  lateral  vibration  is  explored, 
the  author  suggests  that  the  correct  use  of  active  control  will  be  to  only  control  the  lateral 
rocking  modes  caused  by  the  softening  of  the  payload  interface. 
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Appendix  A:  Model  Block  Design  Documentation 

Appendix  A  documents  the  design  of  the  various  Simulink  model  blocks.  The 
emphasis  here  is  on  the  design  for  compatibility  throughout  the  various  models  that  were 
built.  The  author  is  quite  pleased  that  the  final  block  designs  are  fully  compatible  across 
the  several  models  that  were  ultimately  built.  The  concept  used  was  to  assure 
compatibility  so  that  the  researcher  could  quickly  assemble  a  new  model  focused  on 
some  aspect  of  LY  vibration  and  produce  results  that  are  comparable  to  previous  results. 

Focused  model  construction  allows  for  simulation  speed.  A  great  deal  of  effort 
went  into  model  optimization.  Certainly,  there  are  many  places  where  the  models  could 
be  simplified  and  further  enhance  processing  speed.  Potential  areas  for  future  speed 
enhancements  are  converting  the  models  to  State  Flow,  and/or  Real  Time  Workshop. 
Real  Time  Workshop  will  probably  have  the  greatest  impact  on  speed  because  it  converts 
the  Simulink  model  into  executable  code.  Another  area  of  improvement  may  be  to  create 
a  more  detailed  model  of  only  first  stage  flight.  Much  of  the  vibration  loading  on  a 
satellite  occurs  during  first  stage  flight,  and  the  benefits  of  higher  fidelity  could  be 
realized  in  future  controller  work. 

Simulink  Primer 

Three  key  Simulink  features  were  utilized  for  model  construction.  They  are  1) 
sub-systems  and  masks,  2)  vectors,  and  3)  the  workspace.  The  sub-systems  and  masks 
greatly  simplify  model  appearance.  Masks  take  the  sub-system  to  a  greater  level  of 
independence  because  they  allow  the  creation  of  dialog  boxes  from  which  it  is  possible  to 
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enter  variable  values.  This  feature  allows  the  user  to  quickly  change  a  value  of  interest 
and  re-run  a  simulation  with  minimum  fuss. 

The  second  feature  of  note  is  vectors.  Use  of  vectors  requires  dedicated 
documentation  of  signal  lines  or  outputs/inputs.  Vectors  allow  for  greater  use  of  matrix 
functions  and  function  blocks.  The  function  block  reduces  overall  block  and  signal  line 
count  by  combining  several  functions  simultaneously.  Figure  34  emphasizes  this  concept. 
Output  1  and  2  will  result  in  the  same  value,  but  Output  1  required  only  two  blocks  to 
produce.  The  ‘Mux’  block  vectorizes  several  signal  lines.  If  the  Wide  Vector  Lines 
formatting  feature  is  selected,  vectorized  signal  lines  will  appear  as  wider  lines  (see 
Figure  34).  A  ‘Demux’  block  will  de-vectorize  a  signal  line. 


Math 

Figure  34  Example  vector  and  function  model 

Finally,  the  workspace  allows  the  user  to  specify  the  results  of  complex 

calculations  as  variables  to  be  used  by  the  Simulink  model.  This  is  achieved  by  writing 
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and  executing  a  Matlab  ‘.m-file’  or  defining  a  variable  at  the  Matlab  prompt. 
Subsequently  calling  that  variable  from  the  Simulink  model  will  use  its  value  in  the 
model.  This  greatly  simplifies  block  definitions.  A  good  example  is  the  use  of  a  filter  in 
the  noise  model.  Figure  35  is  a  sample  model  of  a  noise  filter.  The  filter  shapes  the 
frequency  content  of  the  band-limited  white  noise  block.  The  filter  numerator  (bp)  and 
denominator  (ap)  are  defined  by  running  a  short  .m-file  prior  to  executing  the  simulation. 
In  this  example,  the  filter  values  are  passed  from  the  workspace  through  a  masked  sub¬ 
system.  Therefore,  the  variable  names  are  called  in  the  mask  dialog  box. 


White  Noise 


Figure  35  Exploded  masked  sub-system  demonstrating  use  of  workspace 

variables 

Excellent  on-line  documentation  is  available  on  the  AFIT  Alpha  computers  at 
file:/local/matlab5/help/helpdesk.html.  The  Helpdesk  has  detailed  summaries  and 
examples  of  Matlab  and  Simulink  functions  or  blocks.  It  also  contains  Portable 
Document  Format  versions  of  the  entire  Simulink  Users  Guide  and  other  Matlab  users 
guides.  Having  a  web  browser  hot-linked  to  this  page  was  essential. 

In  this  discussion  of  the  simulations  and  model  systems  the  block  names  will  be 
italicized. 
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Top-Level  Model 


is  of  the  top-level  model.  With  little  exception,  the  blocks  at  this  level  are  sub¬ 
systems  or  masked  blocks.  This  top-level  view  emphasizes  that  the  clock  drives  the 
model.  The  clock  drives  many  events  of  the  model.  Staging,  gusts  and  other  flight  events 
are  triggered  by  the  clock.  The  model  divides  roughly  into  two  parts  -  environment 
generation  and  dynamic  output.  The  Noise  Generator,  LV  CG  Dynamics,  and  Mass 
Properties  blocks  generate  the  environment  that  excites  the  model  vibration  dynamics. 
The  Axial  or  Lateral  Forcing  Functions  and  the  Axial  or  Lateral  Vibration  Dynamics 
calculate  the  dynamic  output  of  the  model.  The  end  results  are  the  positions,  velocities 
and  accelerations  of  the  LV  components.  The  model  demonstrates  this  flow  from  clock  to 
dynamic  output.  The  blocks  of  the  model  can  be  combined  to  create  new  models. 


Figure  36  Top  level  of  the  Five  DoF  passive  isolation  model 

Environment  generation  in  any  simulation  is  the  same.  The  three  blocks,  LV  CG 

Dynamics,  Mass  Properties  and  the  Noise  Generator,  are  used  universally.  The  blocks 
represent  the  characteristics  of  the  launch  vehicle  and  its  mission  profile.  These  are  the 
only  three  time  dependent  blocks  and  are  the  source  of  the  full  simulation’s  non-linearity. 
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These  will  be  discussed  later  in  this  Appendix.  The  LV  CG  Dynamics  block  is  based  on 
the  gravity-turn  equations  derived  in  Section  2.1.1. 


The  most  basic  model  is  the  five-DoF  simulation  in  Figure  36.  Adding  active 
controls  results  in  Figure  37.  The  LV  model  can  be  radically  simplified  to  use  a  state- 
space  linear  time  invariant  representation  for  rapid  controller  prototyping  as  in  Figure  1 1 . 
Note  that  in  the  state-space  model  the  dynamic  output  does  not  rely  on  the  different 
component  masses  due  to  time-invariance.  This  model  is  only  accurate  at  27  seconds 
after  launch.  The  state-space  block  calls  on  workspace  variables  defined  in  sscont3ppf.m 
(see  Appendix  B)  executed  prior  to  the  simulation. 
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Figure  37  Top  level  of  five  DoF  PID  model 


Figure  38  Top  level  of  the  state  space  estimate  of  the  full  model 
Mass  Properties 

The  very  defining  characteristics  of  a  LV  are  its  masses  and  dimensions.  The 
Mass  Properties  sub-system  mask  defines  these  characteristics.  This  block  is  marked  by 
extensive  use  of  timing  switches  or  step  functions.  A  diagram  is  presented  in  Figure  39. 


Figure  39  Axial  mass  properties  sub-system  diagram 

The  only  input  to  this  block  is  the  simulation  time.  The  time  is  connected  the 

cascade  of  timing  switches  at  the  top  of  the  system  figure.  This  cascade  switches  the  mass 
flow  rate  constants  ( mdot’s )  at  the  staging  events.  The  cascade  is  essentially  an  extended 
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if-then  block  that  determines  the  appropriate  mass  flow  rate  of  the  solid  motor  at  the 
indicated  simulation  time.  Table  3  is  a  listing  of  the  constant  mass  properties  used  to 
define  the  model.  It  should  be  noted  that  these  values  could  be  defined  by  an  .m-file  run 
prior  to  executing  the  simulation.  A  new  set  of  these  values  would  define  a  new  LV. 


Characteristic/Stage 

First 

Second 

Third 

Third/PFJ 

Fourth 

Inert  CM  (m) 

15.6100 

10.8831 

3.9105 

3.2820 

1.2312 

Payload  CM  (m) 

22.8300 

13.8300 

4.9300 

4.9300 

2.2300 

Const  Inert  Inertia  (kg  m2) 

7892606.6 

280826.4 

38742.2 

34039.5 

10.6 

Operating  Stage  Length  (m) 

9 

8.9 

2.7 

2.7 

1.23 

Nozzle  Mass  (kg) 

795.3 

235.4 

118.7 

118.7 

52.8 

Propellant  Density  (kg  m'3) 

1354.0 

1277.1 

960.7 

960.7 

614.0 

Table  3  The  Const  Mass  Props  stored  values 
The  use  of  step  inputs  was  mentioned  previously.  These  are  used  in  lieu  of  timing 

switches  to  specify  certain  mass  change  events  more  explicitly.  The  Simulink  step  input 

function  refers  internally  to  the  simulation  time.  This  creates  some  flexibility  in  defining 

timed  events.  A  switch  is  used  to  define  a  system  state  change.  While  the  step  input 

functions  refer  to  mass  change  events.  System  state  changes  involve  values  that  are 

integrated  and  mass  changes  are  instantaneous  mass  losses. 

Utilizing  the  previous  definitions,  the  cascade  timing  switches  at  the  top  of  the 
diagram  is  more  clearly  understood.  The  empty  stage  masses  and  fairing  mass  appear  as  a 
cascade  of  step  inputs.  These  will  drop  off  at  the  appropriate  moment.  At  the  same  time 
the  switches  determine  which  of  the  vector  of  mdot’s  is  to  be  integrated  by  the  Mass  Flow 
integrator.  These  values  summed  with  the  Payload  Mass  is  the  total  mass  of  the  LY  sent 
to  output  1  -  Tot  Mass.  The  total  mass  is  used  to  calculate  the  mass  values  of  the  Lower 
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and  Upper  model  components  using  the  other  known  system  values.  The  known  values 
are  held  within  the  various  step  functions.  These  include  Pos  Nozzle,  the  Inert  Mass  of 
stage  one,  two,  three  and  the  fairing.  The  payload  mass,  plm,  is  defined  externally  by  the 
user  in  a  dialog  box  that  is  accessed  by  double  clicking  on  the  Mass  Properties  block  at 
the  top  level  of  the  model.  The  mass  of  the  Sub-component,  massl,  is  also  updated  in  the 
same  dialog  box.  It  should  be  obvious  by  inspection  of  the  model  that  the  Lower  and 
Upper  component  mass  are  half  of  the  total  stage  mass  minus  the  nozzle  mass. 

LV  CG  Dynamics 

The  LV  CG  Dynamics  block  creates  the  basis  of  the  acceleration  environment. 
The  block  calculates  Thrust,  Drag,  local  air  density,  pressure,  gravitational  acceleration 
and  the  fundamental  states  of  the  LV’s  center  of  mass.  The  equations  being  modeled  are 
presented  in  Section  2.1.1.  Figure  is  the  block  diagram  of  LV  CG  Dynamics. 

The  diagram  looks  somewhat  complex,  but  it  is  actually  straightforward.  The  four 
integrators  are  the  basis  of  the  block.  These  integrate  the  velocity  ( CG  Velocity ),  altitude 
(CG  Altitude),  range  (CG  Range)  and  flight  path  angle  (CG  Alpha).  CG  Vel  Dot  and 
Alpha  Dot  are  function  blocks.  You’ll  recall  that  the  equations  of  velocity  and  flight  path 
angle  are  somewhat  lengthy.  This  complexity  is  handled  in  the  function  block.  The 
equations  in  the  function  blocks  are  equivalent  to  Equations  (5c)  and  (5d).  The  rest  of  the 
sub-system  brings  together  the  various  environmental  factors  that  make  the  gravity  turn 
equations  very  non-linear. 
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Figure  40  LV  CG  Dynamics  Block  Diagram 

The  Thrust  Program  block  serves  two  purposes.  First,  it  dictates  by  means  of 

switches  the  level  of  thrust  that  the  flight  conditions  dictate.  This  varies  near  sea-level 
because  of  atmospheric  pressure.  Second,  it  calculates  an  exponential  thrust  tail-off 
beginning  seven  seconds  before  burn-out.  Figure  41  is  of  the  Thrust  Program.  There  are 
five  thrust  states.  Four  of  these  states  are  the  stages  and  the  fifth  is  motor  off  or  the  coast 
state. 

Stage  one  (StO  Thrust  in  the  diagram)  and  two  (Stl  Thrust)  thrust  sub-systems  are 
enabled.  This  means  that  the  blocks  are  only  on  when  they  are  operating.  Otherwise  they 
are  turned  off  to  reduce  computational  overhead.  (Note:  Simulink  calculates  the  state  of 
all  blocks  in  a  model  unless  they  are  triggered,  or  enabled.) 
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Figure  41  The  Thrust  Program  Block  Diagram 

The  thrust  tail-off  before  burn-out  is  models  an  effect  all  solid  rocket  motors 

produce.  The  math  model  is  in  Section  2.1.1,  Equation  6. 

Pitch  Program  is  a  simple  sub-system  that  indicates  the  time  and  pitch  angle  for 
the  flight  path  angle  computation.  This  is  similar  to  commanding  thrust  vector  maneuvers 
just  after  lift-off  to  begin  the  gravity  turn  maneuver.  You’ll  note  the  zero  constant  block 
just  below  the  Alpha  Dot  block.  This  is  the  initial  nozzle  deflection  angle.  This  is  fed  to  a 
switch.  The  switch  makes  it  possible  to  keep  Simulink  from  integrating  the  full  flight 
path  angle  equation  until  two  seconds  after  lift-off,  because  this  value  is  undefined  for 
zero  velocity  (See  Equation  5d).  At  which  point  the  notional  rocket  has  cleared  the  tower 
and  can  begin  the  gravity  turn.  The  Pitch  Program  then  commands  a  three  second  ramp 
up  then  a  three  second  ramp  down  to  zero.  The  pitch  program  gain  can  be  set  in  the  LV 
CG  Dynamics  dialog  box. 
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The  Atmosphere  Model  computes  atmospheric  data  for  an  isothermal  atmosphere. 
The  249K  temperature  of  the  model  atmosphere  was  determined  by  integrating  the  layers 
of  the  1976  Standard  Atmosphere  and  averaging  the  temperatures.  Pressure  and  Density 
are  computed  as  appropriate. 

The  user  can  specify  the  ignition  times  of  Third  and  Fourth  stage  as  well  as  the 
pitch  program  gain  and  initial  conditions  for  the  integrators  in  the  dialog  box  for  the  LV 
CG  Dynamics  block.  Specifying  correct  initial  conditions  allows  the  user  to  examine  a 
different  part  of  the  mission  profile  without  having  to  start  at  zero.  Use  the  lvcgdyn.mdl 
model  to  determine  correct  initial  conditions  for  a  given  start  time. 

Inputs  to  the  LV  CG  Dynamics  are  Time,  Sys  Mass,  and  Vibes.  Vibes  originates  at 
the  Noise  Generator  block.  Sys  Mass  originates  at  the  Mass  Properties  block.  Time 
originates  at  the  Clock.  Outputs  are  Thrust,  CG  Vel,  Density  and  Drag  Force. 

Noise  Generator 

The  user  can  specify  the  aeroacoustic  and  motor  noise  power,  magnitude  of 
resonant  burn,  and  frequency  of  resonant  bum  in  the  dialog  box  of  the  Noise  Generator. 
This  block  creates  the  vibration  forces  that  are  the  focus  of  this  thesis.  Figure  42  is  the 
block  diagram  of  the  noise  generator. 
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Figure  42  The  Noise  Generator  Block 

Resonant  bum  is  the  harmonic  oscillation  of  the  burning  gases  inside  the 

combustion  chamber  of  a  solid  rocket  motor  (SRM).  This  problem  is  typical  for  large 
SRMs.  The  on  or  off  state  of  the  resonant  burn  condition  is  modeled  as  a  step  function. 
The  step  function  multiplies  the  sine  function  to  generate  the  oscillating  force.  The 
resonant  bum  condition  only  exists  for  first  stage. 

A  band-limited  white  noise  block  creates  the  random  noise  that  is  the  basis  of  the 
vibration  signal.  The  block  simultaneously  computes  two  values,  one  for  motor  noise  the 
other  for  aeroacoustic  noise.  The  values  are  filtered  to  correct  the  frequency  content 
according  to  the  noisefilt.m  program,  noisefilt.m  must  be  run  at  least  once  prior  to  a 
simulation  session  so  that  the  filter  transfer  function  will  be  in  the  Matlab  workspace. 
These  vibration  signals  are  then  passed  to  the  Combustion  Noise  block  or  multiplied  by 
the  Transonic  Scale  Factor. 
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The  Combustion  Noise  block  amplifies  the  inputted  vibration  signal  by  the  ratio  of 
the  current  thrust  value  to  the  maximum  thrust  of  the  first  stage  motor.  In  this  way  the 
thrust  level  tailors  the  vibration  magnitude.  The  noise  out  is  multiplied  by  the  motor  noise 
power  gain  specified  in  the  dialog  box  for  the  Noise  Generator  block. 


Figure  43  The  Transonic  Scale  Factor  Block 

Referring  to  Section  2.2.2,  Equation  19  we  can  understand  Figure  43.  This  block 

creates  the  profile  of  the  aeroacoustic  vibration  experienced  during  transonic  and  max-q 
flight.  From  the  velocity  input  the  Mach  number  is  calculated,  and  the  dynamic  pressure 
is  calculated  from  the  air  density  input.  The  dynamic  pressure  is  normalized  to 
atmospheres.  A  Mach  effect  value  is  computed,  as  well  as  a  scaled  dynamic  pressure 
effect.  The  sum  of  these  effects  is  multiplied  by  the  atmospheric  density.  The  density  is 
the  final  scale  factor  because  the  effects  of  atmosphere  should  lessen  as  the  air  thins.  The 
resulting  scale  is  used  as  a  gain  on  the  aeroacoustic  vibration  force.  The  vibration  forces 
are  combined  into  a  single  signal  and  passed  to  the  LV  CG  Dynamics  block  and  the  Axial 
Forcing  Functions  block. 
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Axial  Forcing  Functions 


This  relatively  simple  block  creates  the  “right-hand-sides”  of  the  equations  of 
motion.  All  of  the  relevant  forces  are  multiplexed  and  fed  into  one  of  four  force  function 
blocks.  Only  four  are  needed  because  upon  inspection  of  Equations  (10b)  (Upper 
component)  and  (lOe)  (Sub-component)  you’ll  agree  that  the  right  hand  sides  of  these 
equations  are  the  same.  The  forces  input  here  are  Thrust,  Noise,  and  Drag.  The  inverse 
mass  matrix  is  also  input  to  this  block.  The  output  is  formed  into  a  single  multiplex  line 
with  all  five  Axial  Forces.  Figure  44  is  the  block  diagram. 


Figure  44  The  Axial  Forcing  Functions  Block 


Vibration  Dynamics 

This  block  is  ostensibly  the  heart  of  the  simulation.  The  general  appearance  of  this 
block  is  the  classic  state-space  wiring  diagram  (Figure  45).  The  general  layout  should 
need  no  explanation,  but  the  Damping  and  Stiffness  blocks  may  require  discussion. 
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Figure  45  The  Vibration  Dynamics  Block 


Figure  46  The  Damping  Forces  Block 
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Figure  46  is  the  block  diagram  of  Damping.  Damping  and  Stiffness  are  structured 
identically  because  they  are  an  elaboration  of  the  matrix  form  of  the  equation  of  motion. 
Simulink  cannot  pass  matrices,  only  vectors.  This  results  in  some  difficulties  in 
implementing  a  time  varying  mass  matrix.  The  inverse  of  the  mass  values  are  inputted. 
The  total  mass  of  the  rocket  is  separated  from  the  inverse  mass  vector.  The  mass  inverses 
are  then  multiplied  by  a  damping  or  stiffness  vector  appropriate  to  the  row  of  the 
equations  of  motion  the  mass  inverse  belongs.  These  squares  of  the  modal  frequencies  or 
damping  ratios  are  then  multiplied  by  the  position  or  velocity  state  as  appropriate.  These 
are  summed  and  passed  on  as  a  vector  of  damping  or  stiffness  forces. 

PID  Controller 


Gain2  Derivative 


Figure  47  The  PID  Controller  Block 

The  block  diagram  above  is  intuitively  obvious  when  one  considers  Equation  20. 

The  block  is  presented  for  completeness.  The  PPF  controller  is  implemented  using  a 
state-space  block,  provided  by  Simulink. 
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Appendix  B:  .m-file  Analysis  Source  Code 


The  Matlab  source  code  of  the  various  .m-files  developed  for  this  thesis  are 
presented.  The  programs  listed  are  sscont3ppf.m,  freqflux2.m,  rmsg.m,  and  noisefilt.m. 
These  programs  are  fundamental  to  the  analysis  or  production  of  data  by  the  Simulink 
models  used  for  this  study. 


sscont3ppf.m 


%  sscont3ppf  creates  a  state-space  model,  sys,  of  the  launch1 vehicle 
satellite 

%  combination  at  27  seconds.  This  time  is  chosen  to  corespond  to 
transonic 

%  flight.  The  .m-file  then  builds  a  state-space  PPF  controller,  cntrl . 
%  sys  and  cntrl  are  combined  in  closed  loop  form  to  create  sysfeed 
% 

%  Specify  frequency  of  passive  isolator,  f,  and  controller  gain  ,G, 

%  controller  damping,  zc,  and  controller  frequency,  f c . 

% 


%  Declare  PPF  controller  frequencies 
% 

fc=49 . 5 ; 
zc=0 . 05 ; 
f c2=13 . 1 ; 
zc2=0 . 08 ; 
f c3=29 . 0 ; 
zc3=0.05; 

G=  0.002; 

%  System  at  27  sec 

M=  [1912  0  0  0  0  0 ; 0  18330  0  0  0;0  0  18580  0  0;0  0  0  1779  0;  0  0  0  0  1]; 
k=le8 ; 

ki=(2*pi*f) A2*1779; 

K=[4*k  -4*k  0  0  0 ; - 4 * k  5*k  -k  0  0;0  -k  k+ki  -ki  0;0  0  -ki  ki+98696  - 
9 8 6 9 6 ; 0  0  0  -98696  98696] ; 

C= [ 9e4  -9e4  0  0  0;-9e4  9e4*1.5  -45e3  0  0;0  -45e3  54839  -9838.7  0;0  0  - 
9838.7  9844.9  -6.2832;0  0  0  -6.2832  6.2832]; 

F= [ ( 1-M  (1,1) / 57810)  ( 1-M (1,1) /57810)  -M(l, 1) /57810  M(l,l) / 57810 

0; 

-M (2,2) /57810  -M (2,2) /57810  -M (2,2) /57810  M(2,2) / 57 810 

0; 

-M (3,3) /57810  -M(3 , 3 ) /57810  - ( 1-M (3 , 3 ) /57810 )  (1- 

M(3 , 3) /57810)  -G* (2  *pi*f ) A2 ; 
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-1779/57810  -1779/57810  -1779/57810  1779/57810 

G* (2  *pi*f ) A2 ; 

-1/57810  -1/57810  -1/57810  1/57810 

0]  ; 

%  Reduce  System 

T= [ 1  0  0  0; -19120/18330  -18580/18330  -1779/18330  -1/18330;0  1  0  0;0  0  1 
0 ;  0  0  0  1]  ; 

Mr=T ' *M*T; 

Kr=T ' *K*T; 

Cr=T 1  *  C  *  T ; 

Br=T '  *F; 

%create  state  model 
A=zeros (8,8) ; 

Ol=zeros (4,8) ; 

A  ( 1 : 4 , 5 : 8 ) =eye ( 4 ) ; 

A (5 : 8 , 1 : 4) =-inv (Mr) *Kr ; 

01 ( : , 1 : 4) =-inv(Mr) *Kr; 

A (5 : 8 , 5 : 8 ) =-inv(Mr) *Cr ; 

01 (: ,5:8) =-inv (Mr) *Cr; 

B=zeros (8,5) ; 

B  ( 5 : 8 ,  : )=Br; 

C=eye(8, 8) ; 

D=zeros (8,5) ; 
sys=ss (A, B , C , D) ; 

%Create  Controler  Model 


Ac= [ 0  00100 
0  0  0  0  1  0 
0  0  0  0  0  1 

- (2*pi*fc) A2  0  0 

0; 

0  - (2 *pi*fc2 ) ^2  0 

2  *zc2  * (2*pi*fc2 )  0; 

0  0  - (2*pi*f c3 ) A2  0 

-2*zc3* (2*pi*fc3 ) ] ; 

Bc= [ 0  0  0  - ( fc*2*pi) A2  - ( fc2*2*pi) ^2  - ( fc3 *2 *pi) A2 ] 1 ; 
Cc= [1  11000]; 

Dc=0  ; 

cntrl=ss (Ac , Be , Cc , Dc) ; 

sysfeed=feedback (sys , cntrl , [5] , [3] ,+l) ; 


-2*zc* (2*pi*f c)  0 
0 


freqflux2.m 


%  freqflux2  calculates  the  time  varying  modes  of  the  LV  and  satellite 
system 

%  for  a  five  degree  of  freedom  model 
% 

%  Declarations 
% 

1= [ 816 ; 733 ; 808 ; 691] ;  %  length  of  event  counts 

mdo  t=  [-598. 8834  ,*-166.0109; -  37. 48451  ,*-11.3333];  %  mass  flow  rates 
MLi=[26907.7;7127.35;1744.35;518.9] ;  %  Lower  component  masses 
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MUi= [26112 . 35 ; 6891 . 95 ; 1625 . 65 ; 466 . 1] ;  %  Upper  component  masses 
MIi= [18580 ; 4560 ; 1190 ; 985] ;  %  Inert  component  masses 
MPi=1779 ; 

Mli=l; 

kl=4*le8;  %  inter-component  stiffness  KLU 
k2=le8 ;  %  inter-component  stiffness  KUI 

k3= (29*2*pi) A2*MPi;  %  inter-component  stiffness  KIP  in  form 
( f reguencyA2 ) *mass 

k4= (50*2*pi) A2*Mli;  %  inter-component  stiffness  KP1  in  form 
( frequency^  2 ) *mass 
% 

for  x=l:l  %  1  to  number  of  stages 
n=l;  %  event  counts 

t=0 : 0 . 1 : (1 (x) /10- . 1) ;  %  time  vector 

lam=zeros (1 (x) , 4) ;  %  initialize  eigenvalues 

while  n  <=  1 (x) 

M= [MLi (x) +t (n) *mdot (x) /2 ;MUi (x) +t (n) *mdot (x) /2 ] ;  %  calculate 

current  masses 

ML=M(1);  MU=M ( 2 ) ;  MI=MIi (x) ;  MP=MPi;  Ml=Mli; 

%  create  transpostion  matrix 

T= [ 1  0  0  0; -ML/MU  -MI/MU  -MP/MU  -M1/MU;0  1  0  0 ; 0  0  1  0 ; 0  0  0  1 ] ; 
m= [ML  0  0  0  0;0  MU  0  0  0;0  0  MI  0  0;0  0  0  MP  0;0  0  0  0  Ml];  %  mass 
matrix 

%  stiffness  matrix 

K= [kl  -kl  0  0  0 ; -kl  kl+k2  -k2  0  0;0  -k2  k2+k3  -k3  0;0  0  -k3  k3+k4  - 
k4 ; 0  0  0  -k4  k4] ; 

[vec, val] =eig(inv(T' *m*T) *T' *K*T) ;  %  instantaneous 
eigenvalues /vectors 

lam(n,  :  )  =  (  (  (val)  A0 . 5)  *  [1  ;.l;  1 ;  1]  )  1  ;  %  compute  frequencies 
n=n+l;  %  increment  event  counter 
end;  %  while 
% 

%  Plot  output 

figure; plot (t , lam( : , 1 : 4) / (2*pi) ) ; 
title ( 1  Frequency  vs  Time ' ) ; 
xlabel ( 1  Time  ( sec )  ’  )  ; 
ylabel ( 1  Frequency  (Hz)  ' ) ; 
end;  %  for 


rmsg.m 


function  boo=rmsg (n# c) ; 

%computes  rms  in  g's  of  channel,  c,  of  matrix,  n. 
boo=norm(n ( : , c) ) /sqrt (size (n, 1) ) ; 

%end  rmsg 


noisefilt.m 


%  noisefilt  is  to  be  run  prior  to  any  simulation 

%  This  file  shapes  the  frequency  content  of  the  input 

%  band  limited  white  noise  in  the  Noise  Generator  block 

[  z ,  p ,  k] =buttap ( 4 ) ;  %  Build  fourth  order  Butterworth  Prototype  Filter 

[b, a] =zp2tf (z,p, k) ;  %  Convert  to  transfer  function  format 

[bp, ap] =lp2hp (b, a, 20*pi ) ;  %  Convert  filter  from  lowpass  to  high  pass 

form 

clear  a;  clear  b;  clear  k;  clear  z;  clear  p;  %  clear  extraneous 
variables 


Appendix  C:  Model  Source  Code 


The  Simulink  source  code  of  the  PID  hybrid  model  is  provided. 


Model  { 


Name 

"  fi  vedofhybridPIDFL' 

Version 

2.09 

SimParamPage 

Diagnostics 

SampleTimeCoIors 

off 

InvariantConstants 

off 

WideVectorLines 

on 

ShowLineWidths 

off 

PaperOrientation 

landscape 

PaperType 

usletter 

PaperUnits 

inches 

StartTime 

"OX)" 

StopTime 

"90" 

Solver 

ode45 

RelTol 

"le-3" 

AbsTol 

"le-6" 

Refine 

Mr 

MaxStep 

"auto" 

InitialStep 

"auto" 

FixedStep 

"1/800" 

MaxOrder 

5 

OutputOption 

RefineOutputTimes 

OutputTimes 

"[  1" 

LoadExtemallnput 

off 

Extemallnput 

"ft  u]" 

SaveTime 

off 

TimeSaveName 

"tout" 

SaveState 

off 

StateSaveName 

"xout" 

SaveOutput 

off 

OutputSaveName 

"yout” 

LoadlnitialState 

off 

InitialState 

"xlnitial" 

SaveFinalState 

off 

FinalStateName 

"xFinal" 

LimitMaxRows 

off 

MaxRows 

"1000" 

Decimation 

"1" 

AlgebraicLoopMsg 

warning 

MinStepSizeMsg 

warning 

UnconnectedlnputMsg 

none 

UnconnectedOutputMsg 

none 

UnconnectedLineMsg 

warning 

Consistencychecking 

off 

ZeroCross 

on 

SimulationMode 

normal 

RTWSystemTargetFile 

"grt.tlc” 

RTWInlineParameters 

off 

RTWRetainRTWFile 

off 

RTWTemplateMakefile 

"grt_unix.tmf 

RTWM  akeCommand 

"make_rtw" 

RTWGenerateCodeOnly 

off 

ExtModeMexFile 

"ext_comm" 

ExtModeBatchMode 

off 
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BlockDefauIts  { 


Orientation 

right 

ForegroundColor 

black 

BackgroundColor 

white 

DropShadow 

off 

NamePlacement 

normal 

FontName 

"Helvetica" 

FontSize 

10 

FontWeight 

normal 

FontAngle 

normal 

ShowName 

on 

} 

AnnotationDefaults  { 


HorizontalAlignment 

center 

VerticalAlignment 

middle 

ForegroundColor 

black 

BackgroundColor 

white 

DropShadow 

off 

FontName 

"Helvetica" 

FontSize 

10 

FontWeight 

normal 

FontAngle 

i 

normal 

I 

LineDefaults  { 

FontName 

"Helvetica" 

FontSize 

9 

FontWeight 

normal 

FontAngle 

\ 

normal 

; 

System  { 

Name 

"fivedofhybridPIDFL" 

Location 

[227,588,  1257,935] 

Open 

on 

ScreenColor 

white 

Block  { 

BlockType 

Subsystem 

Name 

"Axial  Forcing\nFunctions" 

Ports 

[5, 1,0,0,  0] 

Position 

[615,  27,715,  223] 

NamePlacement 

alternate 

ShowPortLabels 

on 

System  { 

Name 

"Axial  Forcing\nFunctions' 

Location 

[141,34,707,317] 

Open 

off 

ScreenColor 

white 

Block  { 

BlockType 

Inport 

Name 

"Thrust" 

Position 

[15,  18,  45,  32] 

Port 

"1" 

PortWidth 

M" 

SampleTime 

i 

”_r 

/ 

Block  { 

BlockType 

Inport 

Name 

"Noise" 

Position 

[15,  50,  45, 70] 

Port 

"2" 

PortWidth 

"-1" 

SampleTime 

i 

"-r 

/ 

Block  { 

BlockType 

Inport 

Name 

"Drag" 

Position 

[15,88,  45,102] 

Port 

"3" 

PortWidth 

"_r 
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SampleTime 

} 

Block  { 
BlockType 
Name 
Position 
Port 

PortWidth 

SampleTime 

> 

Block  { 
BlockType 
Name 
Position 
Port 

PortWidth 

SampleTime 

} 

Block  { 
BlockType 
Name 
Ports 
Position 
Outputs 

} 

Block  { 
BlockType 
Name 
Ports 
Position 
Inputs 

} 

Block  { 
BlockType 
Name 
Ports 
Position 
Inputs 

} 

Block  { 
BlockType 
Name 
Position 
Expr 


"-1" 


Inport 
"Mass" 
[15,  150,35,  170] 
"4" 

"-1" 


Inport 
"Isolator" 
[120,  235, 140,  255] 
"5" 

M" 

"-l" 


Demux 
"Demux" 
[1,6,  0,  0,  0] 
[50,  117,  90, 198] 

"6" 


Mux 

"Force  Inputs" 
[8,  1, 0,  0,  0] 
[210, 79,  250,  161] 

"8" 


Mux 

"Forcing\nVector\n  Inputs" 
[5, 1,0,  0,0] 

[435,91,470,  149] 

"5" 


Fen 

"I  Force" 

[305, 105,  365, 135] 

"-u[7]*(u[l]+u[2])+(u[7]-u[5])*(u[4]+u[3])-u[" 

”9]/u[5]" 


Block  { 
BlockType 
Name 
Position 
Expr 

} 

Block  { 
BlockType 
Name 
Position 
Expr 

} 

Block  { 
BlockType 
Name 
Position 
ShowName 

} 

Block  { 
BlockType 
Name 
Position 


Fen 

"L  Force" 

[305,15,  365,45] 

"(u[6]-u[7])*(u[l]+u[2])+u[7]*(u[4]+u[3])" 


Fen 

"P  Force" 

[305,  155,  365,  185] 

"-u[7]*(u[l]+u[2]-u[3]-u[4])+u[9]/u[8]" 


Terminator 
"Terminator23" 
[140,  142,  160,  158] 
off 


Terminator 
MTerminator25" 
[140,  188,  160,  202] 
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ShowName 


off 


} 

Block  { 

BlockType 

Name 

Position 

Expr 

} 

Block  { 

BlockType 

Name 

Position 

Port 

OutputWhenDisabled 

InitialOutput 

} 

Line  { 

SrcBlock 

SrcPort 

Points 

DstBlock 

DstPort 


Fen 

"U  Force" 

[305,  60,  365,  90] 

"-u[7]*(u[l]+u[2]-u[4]-u[3])” 


Outport 

"Axial\nForce\nVector" 


[490,  113,520,  127] 
"1" 
held 
"0" 


"P  Force" 

1 

[50,  0] 

"Forcing\nVector\n  Inputs" 
4 


} 

Line  { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 


"Demux" 

5 

[90,  0;  0,  -35] 
"Force  Inputs" 
7 


} 

Line  { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 


"Isolator" 

1 

[45,  0;  0,  -90] 
"Force  Inputs" 
8 


"Forcing\nVector\n  Inputs" 
1 

"Axial\nForce\nVector" 

1 


Line  { 
Labels 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 


[1,0] 

"L  Force" 

1 

[50,  0] 

"Forcing\nVector\nInputs" 

1 


} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 


"I  Force" 

1 

"Forcing\nVector\n  Inputs" 
3 


Line  { 
SrcBlock 
SrcPort 
Points 
Branch  { 
Points 
DstBlock 
DstPort 

} 

Branch  { 
DstBlock 


"U  Force" 

1 

[30, 0;  0,  35] 


[0,  30] 
5 


"Forcing\nVector\nInputs" 


MForcing\nVector\nInputs" 
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DstPort 

} 

) 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 


2 


"Demux” 

6 

"Terminator25" 

1 


"Demux" 

3 

"Terminator23" 

1 


"Demux" 

2 

[15,  0;  0,-10] 
"Force  Inputs" 
5 


"Demux" 

4 

[25,  0;  0,  -50] 
"Force  Inputs" 
4 


} 

Line  { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
Points 
Branch  { 
Points 
DstBlock 
DstPort 


"Demux" 

1 

[0,  10;  20,  0;  0,  5] 
"Force  Inputs" 

6 


"Mass" 

1 

"Demux" 

1 


"Force  Inputs" 
1 

[20,  0] 

[0,50] 


1 


"P  Force" 


Branch  { 
DstBlock 
DstPort 

} 

Branch  { 
Points 
Branch  { 
Points 
DstBlock 
DstPort 

} 

Branch  { 
DstBlock 
DstPort 

} 

} 


"I  Force" 

1 


[0,  -45] 

[0,  -45] 

"L  Force" 
1 


"U  Force" 

1 
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Line  { 
SrcBIock 
SrcPort 
Points 
DstBlock 
DstPort 

} 

Line  { 
Labels 
SrcBIock 
SrcPort 
Points 
DstBlock 
DstPort 

} 

Line  { 
SrcBIock 
SrcPort 
Points 
DstBlock 
DstPort 

} 

} 

} 

Block  { 

BlockType 

Name 

Position 

Location 


"Thrust'’ 

1 

[145,0] 

"Force  Inputs" 

1 


[2,0] 

"Noise" 

1 

[140,  0;  0,  35] 
"Force  Inputs" 
2 


"Drag" 

1 

[40, 0;  0, 10] 
"Force  Inputs" 
3 


Clock 

"Clock" 

[15,  55,  35,  75] 
[21,63,  128,  108] 


Block  { 
BlockType 
Name 
Ports 
Position 
Outputs 

} 

Block  { 
BlockType 
Name 
Ports 
Position 

NamePlacement 

ShowPortLabels 

MaskType 

MaskDescription 


MaskHelp 


MaskPromptString 


MaskStyleString 

MaskVariables 

MasklconFrame 

MasklconOpaque 

MasklconRotate 

MasklconUnits 

MaskValueString 

System  { 

Name 


Demux 

"Demux" 

[1,5,  0,  0,  0] 

[910,  27,  945, 63] 
"5" 


Subsystem 
"LV  CG  Dynamics" 

[3,  4,  0,  0,  0] 

[255,  32,  360,  173] 

alternate 

on 

"LV  CG  Dynamics" 

"Computes  flight  path  of  launch  vehicle  CG. " 

"Vary  final  orbit  altitude  by  adjusting  Payload  " 
"Mass,  Pitch  Program  Gain  and  Fourth  Stage  " 
"Ignition  Time.\n\nNote:  4th  Stage  ignition  " 

"time  must  be  at  least  80.7  sec  after  3rd  " 

"Stage  Ignition." 

"Circular  Orbits  are  possible  with  the  following  " 
"combinations :\n3d  Ig=  166.8  \n4th  Ig=635  " 

"  \nPitch  Gain  =0.133127" 

"2nd  Stage  Ignition  Timel3rd  Stage  ignition  " 
"TimelPitch  GainIRange  Initial  ConditionslAltitud" 
"e  Initial  ConditionIVelocity  Initial  Condition  " 
"IFlight  Path  Angle  Initial  Condition" 

"edit, edit, edit, edit, edit, edit, edit" 

"st2ig=@l;st3ig=@2;pitch=@3;mought=@4;hnought=@5" 

";vnought=@6;anought=@7;" 

on 

on 

none 

autoscale 

"166.8163510. 133 1271010 IOIpi/2" 

"LV  CG  Dynamics" 
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Location 

Open 

ScreenColor 
Block  { 
BlockType 
Name 
Position 
Port 

PortWidth 

SampleTime 

} 

Block  { 
BlockType 
Name 
Position 
Port 

PortWidth 

SampleTime 


[157, 141,928,  651] 
off 
white 

Inport 
"Time” 
[20, 148,  50,  162] 
"1" 

”-1” 

”-1” 


Inport 

"Sys  Mass” 
[20, 113,  50,  127] 

"2" 

n-r 

"-r 


Block  { 
BlockType 
Name 
Position 
Port 

PortWidth 

SampleTime 


Inport 
"Vibes” 
[20,  223,  50,  237] 
"3” 

"-l" 

"-r 


Block  { 
BlockType 
Name 
Position 
NamePIacement 
Expr 


} 

Block  { 

BlockType 

Name 

Ports 

Position 

Orientation 

NamePIacement 

ShowPortLabels 

MaskType 

MaskDescription 


MasklconFrame 
MasklconOpaque 
MasklconRotate 
MasklconUnits 
System  { 

Name 

Location 

Open 

ScreenColor 
Block  { 
BlockType 
Name 
Position 
Port 

PortWidth 

SampleTime 

} 

Block  { 
BlockType 
Name 
Position 


Fen 

"Alpha  Dot" 

[210,320,  270,  350] 
alternate 

"-(u  [  1  ]  *u[4]  -(u[5]  *u  [3]  *u(6)+u  [4]  *u  [3]  A2)/(63 " 
"78135+u[2]))/u[3]" 


Subsystem 
"Atmospere  Model" 

[1,  1,  0,  0,  0] 

[430, 77,  535,  103] 
left 

alternate 

on 

"Atmosphere  Model" 

"Calculates  exponential  pressure  and  density  " 

"based  on  isothermal  atmosphere.  T=249K.\n\nu" 
"( 1  )=pressure\nu(2)=density" 
on 
on 
none 

autoscale 


"Atmospere  Model" 

[647,715,982,847] 

off 


white 


Inport 

"Altitude" 

[25,  58,  55,  72] 


Fen 

"Density" 

[120,  50,  180,  80] 
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Expr 

"1.2250  * 

exp(-(1.3728e-4)  *  u)1 

Block  { 

BlockType 

Outport 

Name 

"pres/dens 

Position 

[260,  58,  290,  72] 

Port 

"r 

OutputWhenDisabled 

held 

InitialOutput 

"0" 

Line  { 

SrcBlock 

"Density" 

SrcPort 

1 

DstBlock 

"pres/dens" 

DstPort 

1 

/ 

Line  { 

Labels 

[1.0] 

SrcBlock 

"Altitude" 

SrcPort 

1 

DstBlock 

"Density" 

DstPort 

1 

} 

} 

} 

Block  { 


BlockType 

Integrator 

Name 

"CG  Alpha" 

Ports 

[1,  1,0,  0,  0] 

Position 

[350,  325,  380,  355] 

ExtemalReset 

none 

InitialConditionSource  internal 

InitialCondition 

"anought" 

LimitOutput 

off 

UpperSaturationLimit 

"inf 

LowerSaturationLimit 

"-inf’ 

ShowSaturationPort 

off 

ShowStatePort 

off 

AbsoluteTolerance 

} 

Block  { 

"auto" 

BlockType 

Integrator 

Name 

"CG  Altitude' 

Ports 

[L  1,0,  0,  0] 

Position 

[595,  310,  625,  340] 

ExtemalReset 

none 

InitialConditionSource  internal 

InitialCondition 

"hnought" 

LimitOutput 

off 

UpperSaturationLimit 

"inf 

LowerSaturationLimit 

"-inf 

ShowSaturationPort 

off 

ShowStatePort 

off 

AbsoluteTolerance 

"auto" 

1 

Block  { 

BlockType 

Integrator 

Name 

"CG  Range" 

Ports 

[1,  1,0,  0,  0] 

Position 

[595,  370,  625,  400] 

ExtemalReset 

none 

InitialConditionSource  internal 

InitialCondition 

"mought" 

LimitOutput 

off 

UpperSaturationLimit 

"inf 

LowerSaturationLimit 

"-inf 

ShowS  aturationPort 

off 

ShowStatePort 

off 

C-8 


AbsoluteTolerance  "auto1' 


} 

Block  { 
BlockType 
Name 
Position 
Expr 


Fen 

"CG  Vel  Dot" 

[350,  155,410,  185] 

"(u[3]+u[10]-u[ll]-u[4]-(u[2]*u[6]*u[7]+(u[l]" 

"*u[9]*u[8]-u[2]*u[7]*u[l]A2)/(6378135+u[5]))" 

')/u[2]" 


} 

Block  { 

BlockType 
Name 
Ports 
Position 
ExtemalReset 
InitialConditionSource  internal 
InitialCondition  "vnought" 
LimitOutput 
UpperSaturationLimit 
LowerSaturationLimit 
ShowSaturationPort  off 
ShowStatePort 
AbsoluteTolerance  "auto" 


Integrator 
"CG  Velocity" 
[1, 1,0,  0,0] 
[445,  155,475,  185] 
none 


off 

"inf 

"-inf 

off 


} 

Block  { 
BlockType 
Name 
Position 
ShowName 
Value 


Constant 
"Constant" 
[225, 370,  245,  390] 
off 
"0" 


} 

Block  { 
BlockType 
Name 
Position 
Orientation 
NamePlacement 
Expr 

} 

Block  { 
BlockType 
Name 
Ports 
Position 
ShowName 
Inputs 

} 

Block  { 
BlockType 
Name 
Ports 
Position 
ShowName 
Inputs 

} 

Block  { 
BlockType 
Name 
Ports 
Position 
Orientation 
ShowName 
Inputs 

} 

Block  { 
BlockType 
Name 


Fen 

"Drag" 

[275,  85,  335,  115] 
left 

alternate 

"0.08*u[l]*u[2]A2" 


Mux 

"Mux" 

[6,  1,  0,  0,  0] 
[160,  308, 195,  362] 
off 

"6" 


Mux 

"Muxl" 

[10,  1,0,  0,  0] 
[300,  125,335,215] 
off 

"10" 


Mux 

"Mux2" 

[2,  1,  0,  0,  0] 
[355,  82,  390,  113] 
left 
off 

"2" 


Subsystem 
"Pitch  Program" 


C-9 


Ports 
Position 

ShowPortLabels 
System  { 

Name 

Location 

Open 

ScreenColor 

Block  { 

BlockType 

Name 

Position 

Port 

PortWidth 

SampleTime 

} 

Block  { 

BlockType 

Name 

Position 

ShowName 

Value 

} 

Block  { 

BlockType 

Name 

Position 

ShowName 

Value 

} 

Block  { 

BlockType 

Name 

Ports 

Position 

ShowPortLabels 

System  { 

Name 

Location 

Open 

ScreenColor 
Block  { 
BlockType 
Name 
Ports 
Position 


[1,  1,0,  0,  0] 

[285,  421,355,469] 
on 

"Pitch  Program" 

[253,  605,  558,  793] 
off 


"Time" 


*r 


white 

Inport 

[25,  53,  55,  67] 
"-1" 


H-r 


"Constant" 


"1” 


Constant 

[30,  90,  50, 110] 
off 


Constant 
"Constant  1" 

[50,  25,  70,  45] 
off 

"0" 


Subsystem 
"Pitch  Program" 

[0,  1, 1,0,  0] 

[185,  119,  225, 161] 
on 

"Pitch  Program" 
[664,  675,1030,  882] 
off 
white 

EnablePort 

"Enable" 

[0,  0, 0,  0,  0] 

[30,  20,  50,  40] 


States WhenEnabling  held 
ShowOutputPort  off 


Block  { 
BlockType 
Name 
Ports 
Position 
ShowName 
Inputs 

} 

Block  { 
BlockType 
Name 
Position 
Gain 

} 

Block  { 
BlockType 
Name 
Ports 
Position 


Sum 

"Event2\nTvectorl " 
[3,  1,  0,  0,  0] 

[160,  102, 195, 138] 
off 

"+++" 


Gain 

"Gain" 

[215,  105,275,  135] 
"pitch" 


Reference 

"Rampl" 

[0, 1,  0,  0,  0] 

[90,  75,  120,  105] 


C-10 


ShowName 

SourceBlock 

SourceType 

slope 

start 

XO 

} 

Block  { 
BlockType 
Name 
Ports 
Position 
ShowName 
SourceBlock 
SourceType 
slope 
start 
XO 


off 

"simulink/Sources/Ramp" 

"Ramp" 

"1/3" 

"2.0" 

"■r 


Reference 

"Ramp2" 

[0,  1,  0,  0,  0] 

[50,  105,  80,135] 
off 

"simulink/Sources/Ramp" 

"Ramp" 

"-2/3" 

"5" 

1 


Line  { 
SrcBIock 
SrcPort 
DstBlock 
DstPort 

} 

Line  { 
SrcBIock 
SrcPort 
DstBlock 
DstPort 

} 

Line  { 
SrcBIock 
SrcPort 
Points 
DstBlock 
DstPort 

} 

Line  { 
SrcBIock 
SrcPort 
DstBlock 
DstPort 


"Gain" 

1 

"Pitch" 

1 


"Event2\nTvectorl" 

1 

"Gain" 

1 


"Rampl" 

1 

[0,  20] 

"Event2\nTvectorl " 

1 


"Ramp2" 

1 

"Event2\nTvectorl " 
2 


} 

Line  { 

SrcBIock  "Ramp3" 

SrcPort  1 

Points  [0,  -20] 


C-ll 


DstBlock 

DstPort 


"Event2\nTvectorl 

3 


} 

} 

} 

Block  { 

BlockType 

Name  "Switch" 

Position 

ShowName 

Threshold 

} 

Block  { 

BlockType 

Name  "Switch  1" 

Position 

ShowName 

Threshold 

} 

Block  { 

BlockType 

Name  "Pitch" 

Position 

Port  "1" 

OutputWhenDisabled  held 
InitialOutput  "0" 

} 

Line  { 

SrcBIock 

SrcPort 

DstBlock 

DstPort 

} 

Line  { 

SrcBIock 

SrcPort 

Points  [20, 0] 

DstBlock 

DstPort 


Switch 

[90, 45,  120,  75] 
off 
"8" 


Switch 

[150,  80,180,110] 
off 
"2" 


Outport 

[250,  133,  280,  147] 


"Pitch  Program" 

1 

"Pitch" 

1 


"Switchl" 

1 

"Pitch  Program" 
enable 


} 

Line  { 

SrcBIock 
SrcPort 
Points 
Branch  { 

Points 

DstBlock 

DstPort 

} 

Branch  { 

DstBlock 

DstPort 


} 


} 

Line  { 

SrcBIock 
SrcPort 
Points 
Branch  { 

Points 

DstBlock 

DstPort 


} 

Branch  { 


DstBlock 

DstPort 


[15,  0] 


"Time" 

1 


[0,  35] 

"Switchl" 

2 


"Switch" 

2 


[0,  0] 


"Constantl" 

1 


[50,  0;  0,  70] 

"Switchl" 

3 


"Switch" 

1 


} 


C-12 


Line  { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 

} 

} 

} 

Block  { 
BlockType 
Name 
Ports 
Position 
ShowName 
Inputs 

} 

Block  { 
BlockType 
Name 
Ports 
Position 
ShowName 
Inputs 

} 

Block  { 
BlockType 
Name 
Ports 
Position 
ShowName 
Inputs 

} 

Block  { 
BlockType 
Name 
Position 
ShowName 
Threshold 

} 

Block  { 
BlockType 
Name 
Position 
ShowName 

} 

Block  { 
BlockType 
Name 
Ports 
Position 

ShowPortLabels 
System  { 

Name 

Location 

Open 

ScreenColor 
Block  { 
BlockType 
Name 
Position 


’’Constant" 

1 

[10,  0;  0,  -30] 

"Switch" 

3 


"Switch" 

1 

[10,  0] 

"Switch  1" 

1 


Product 

"Product" 

[2,  1,0,  0,0] 
[535,  307,  565,  338] 
off 

"2" 


Product 
"Product  1" 
[2,  1,0,  0,0] 
[535,  367,  565,  398] 
off 

"2" 


Sum 

"Sum" 

[2,  1,0,  0,0] 
[410,332, 440,  363] 
off 


Switch 
"Switch" 
[295,  340,  325,  370] 
off 
"2.0" 


Terminator 
"Terminator26" 
[660,  375,  680,  395] 
off 


Subsystem 

"Thrust\nProgram" 

[2, 1,  0,  0,  0] 

[170,  141,220,  199] 
on 

"Thrust\nProgram" 

[144,  599, 746,  891] 
off 

white 

Inport 

"Time" 

[45,  33,  75,  47] 


Port 

PortWidth 

SampleTime 

} 

Block  { 
BlockType 
Name 
Position 

NamePlacement 

Port 

PortWidth 

SampleTime 

} 

Block  { 
BlockType 
Name 
Position 
ShowName 
Value 

} 

Block  { 
BlockType 
Name 
Position 
ShowName 
Value 

} 

Block  { 
BlockType 
Name 
Position 
Value 


T 

"-1" 


M" 


Inport 

"Altitude” 

[25,183,  55, 197] 

alternate 

"2” 

"-1" 

”-r 


Constant 

"Constant" 

[75,215,95,  235] 
off 


Constant 
"Constant  1” 

[75,195,95,215] 

off 


Constant 
"Motor  Off 

[120,  143,  180,  167] 

"0" 


} 

Block  { 

BlockType 

Name 

Ports 

Position 

NamePlacement 
ShowPortLabels 
System  { 

Name 

Location 

Open 

ScreenColor 

Block  { 
BlockType 
Name 
Position 
Port 

PortWidth 

SampleTime 

} 

Block  { 
BlockType 
Name 
Ports 
Position 
Inputs 

} 

Block  { 
BlockType 
Name 
Ports 
Position 
Inputs 


Subsystem 
"St  2  Thrust" 

[1,  1,0,  0,  0] 

[120,  104,  180,  136] 

alternate 

on 

"St  2  Thrust" 
[221,351,758,518] 
off 
white 

Inport 

"Time" 

[85,  15,  105,  35] 

"1" 

"-1" 

M" 


Mux 

"Mux" 

[2,  1,  0,  0,  0] 

[190, 17,  225,53] 

"2" 


Product 

"Product" 

[2,  1,0,0,  0] 

[355,  27,  385,  58] 

"2" 


Block  { 


C-14 


BlockType 

Constant 

Name 

"St2  Thrust" 

Position 

[245,  108,  305,  132] 

Value 

} 

Block  { 

"1 80000" 

BlockType 

Step 

Name 

"Step" 

Position 

[140,  30,  170,  60] 

Time 

"[st2ig+73  st2ig+80]" 

Before 

"[0  D" 

After 

i 

’’[1  0]” 

I 

Block  { 

BlockType 

Switch 

Name 

"Switch" 

Position 

[410,  65,  440,  95] 

Threshold 

} 

Block  { 

"st2ig+73" 

BlockType 

Fen 

Name 

"Thrust  Tail  off 

Position 

[250,  20,  310,  50] 

Expr 

\ 

"u[2]  *exp(-(u[  l]-(st2ig+73)))' 

; 

Block  { 

BlockType 

Outport 

Name 

"St2Thrust" 

Position 

[485,  70,  505,  90] 

Port 

"1" 

OutputWhenDisabled 

held 

InitialOutput 

} 

Line  { 

"0" 

SrcBlock 

"St2  Thrust" 

SrcPort 

1 

Points 

Branch  { 

[30,  0] 

DstBlock 

"Product" 

DstPort 

2 

) 

Branch  { 

Points 

[55,  0] 

DstBlock 

"Switch" 

DstPort 

} 

3 

} 

Line  { 

SrcBlock 

"Switch" 

SrcPort 

1 

DstBlock 

"St2Thrust" 

DstPort 

1 

} 

Line  { 


SrcBlock 

"Mux" 

SrcPort 

1 

DstBlock 

"Thrust  Tail  off 

DstPort 

} 

Line  { 

1 

SrcBlock 

"Step" 

SrcPort 

1 

DstBlock 

"Mux" 

DstPort 

} 

Line  { 

2 

SrcBlock 

"Thrust  Tail  off 

SrcPort 

1 

C-15 


DstBlock 

DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
Points 
Branch  { 
DstBlock 
DstPort 

} 

Branch  { 
Points 
DstBlock 
DstPort 

} 


"Product" 

1 


"Product" 

1 

[5,0] 

"Switch" 

1 


"Time" 

1 

[15,0] 

"Mux" 

1 


[0,55] 

2 


"Switch" 


} 

Block  { 

BlockType  Subsystem 

Name  "St  3  Thrust" 

Ports  [1,1, 0,0,0] 

Position  [120,  54,  180,  86] 

NamePIacement  alternate 


ShowPortLabels 
System  { 

Name 

Location 

Open 

ScreenColor 

Block  { 
BlockType 
Name 
Position 
Port 

PortWidth 

SampleTime 

} 

Block  { 
BlockType 
Name 
Ports 
Position 
Inputs 

} 

Block  { 
BlockType 
Name 
Ports 
Position 
Inputs 

} 

Block  { 
BlockType 
Name 
Position 
Value 

} 

Block  { 
BlockType 


"St  3  Thrust" 
[23.9, 438, 776,  605] 
off 
white 

Inport 

"Time" 

[85,  15,  105,  35] 

"r 

"-i" 


Mux 

"Mux" 

[2,  1, 0,  0,  0] 

[190, 17,  225,  53] 

"2" 


Product 

"Product" 

[2,  1, 0,  0,  0] 

[355, 27,385,  58] 

"2" 


Constant 
"St3  Thrust" 
[250,  108,  310,  132] 
"35000" 


Step 


C-16 


Name 

"Step" 

Position 

[140,  30, 170,  60] 

Time 

"[st3ig+62  st3ig+69]" 

Before 

"[0  1]" 

After 

} 

Block  { 

"[1  03- 

BlockType 

Switch 

Name 

"Switch" 

Position 

[410,  65,  440,  95] 

Threshold 

} 

Block  { 

"st3ig+62" 

BlockType 

Fen 

Name 

"Thrust  Tail  off 

Position 

[250,  20,  310,  50] 

Expr 

} 

Block  { 

"u[2]*exp(-(u[l]-(st3ig+62)))' 

BlockType 

Outport 

Name 

"St3Thrust" 

Position 

[485,  70,  505,  90] 

Port 

"1" 

OutputWhenDisabled 

held 

InitialOutput 

"0" 

; 

Line  { 

SrcBlock 

"St3  Thrust" 

SrcPort 

1 

Points 

Branch  { 

[25,0] 

DstBlock 

"Product" 

DstPort 

} 

Branch  { 

2 

Points 

[55,0] 

DstBlock 

"Switch" 

DstPort 

} 

3 

} 

Line  { 

SrcBlock 

"Switch" 

SrcPort 

1 

DstBlock 

"St3Thrust" 

DstPort 

} 

Line  { 

1 

SrcBlock 

"Mux" 

SrcPort 

1 

DstBlock 

"Thrust  Tail  off 

DstPort 

} 

Line  { 

1 

SrcBlock 

"Step" 

SrcPort 

1 

DstBlock 

"Mux" 

DstPort 

} 

Line  { 

2 

SrcBlock 

"Thrust  Tail  off 

SrcPort 

1 

DstBlock 

"Product" 

DstPort 

} 

Line  { 

1 

SrcBlock 

"Product" 

SrcPort 

1 

Points 

[5,  0] 

C-17 


DstBIock 

DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
Points 
Branch  { 
DstBIock 
DstPort 

} 

Branch  { 
Points 
DstBIock 
DstPort 


"Switch" 

1 


"Time" 

1 

[15,  0] 

"Mux" 

1 


[0,  55] 
2 


"Switch" 


} 

} 

} 

} 

Block  { 

BlockType 

Name 

Ports 

Position 

ShowPortLabels 

System  { 

Name 

Location 

Open 

ScreenCoIor 
Block  { 
BlockType 
Name 
Position 
Port 

PortWidth 

SampleTime 


Subsystem 
"StO  Thrust" 

[2,  1, 1,0,0] 

[210,  226,  340,  259] 
on 


"StO  Thrust" 
[297,547,819,  829] 
off 
white 

Inport 
"Altitude" 
[30,  200,  50, 220] 

"1" 

"_r 

"-i" 


Block  { 
BlockType 
Name 
Position 
Port 

PortWidth 

SampleTime 

} 


Inport 

"Time" 

[30,  80,  50,  100] 
"2" 

"-1" 

"-1" 


Block  { 
BlockType 
Name 
Ports 
Position 


EnablePort 

"Enable" 

[0,  0,  0,  0,  0] 

[30,  25,  50, 45] 


StatesWhenEnabling  held 
ShowOutputPort  off 


Block  { 

BlockType 

Name 

Ports 

Position 

ShowPortLabels 

MaskType 

MaskDescription 


MasklconFrame 

MasklconOpaque 

MasklconRotate 


Subsystem 
"Atmospere  Model  1" 

[1,  1,0,  0,  0] 

[70,  192,  175,  228] 
on 

"Atmosphere  Model" 

"Calculates  exponential  pressure " 

"and  density  based  on  isothermal " 
"atmosphere.  T=249K.\n\nu(  1  )=pressure" 
"\nu(2)=density" 
on 
on 
none 


C-18 


MasklconUnits 
System  { 

Name 

Location 

Open 

ScreenColor 
Block  { 
BlockType 
Name 
Position 
Port 

PortWidth 

SampleTime 


autoscale 


"Atmospere  Model  1” 

[167,351,502,483] 

off 

white 


"Altitude" 

"P 

"-P 


Inport 

[25,  58,  55,  72] 
"-1" 


Block  { 

BlockType 

Name 

Position 

Expr 

} 

Block  { 

BlockType 

Name 

Position 

Port 

OutputWhenDisabled 

InitialOutput 

} 

Line  { 

Labels 

SrcBlock 

SrcPort 

Points 

DstBlock 

DstPort 


Fen 

"Pressure" 

[120,  25,  180,  55] 

"  101325*exp((-l  ,3728e-4)*u)" 


Outport 

"pres" 

[280,58,310,72] 

"P 

held 

"0" 


[1.0] 

"Altitude" 

1 

[15,0;  0,  -25] 

"Pressure" 

1 


} 

Line  { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 

} 

} 

} 

Block  { 
BlockType 
Name 
Ports 
Position 
Inputs 

} 

Block  { 
BlockType 
Name 
Ports 
Position 
Inputs 

} 

Block  { 
BlockType 
Name 
Position 
Expr 

) 

Block  { 
BlockType 
Name 


"Pressure" 

1 

[20,  0;  0,  15;  55,  0;  0,10] 
"pres" 

1 


Mux 

"Mux" 

[2,  1,0,0,  0] 

[135,  82,  170,  118] 

"2" 


Product 

"Product" 

[2,  1,0,0,  0] 

[280,  92,  310,  123] 

"2" 


Fen 

"StO  Thrust" 

[195,  195,  255,  225] 

"1313900+(1 -u/101325)*336100" 


Step 

"Step" 


C-19 


Position 

[85,95,  115,  125] 

Time 

"[74.5  81.53" 

Before 

"[0  1]" 

After 

} 

Block  { 

"[1  0]” 

BlockType 

Switch 

Name 

"Switch” 

Position 

[335,  160,  365,  190] 

Threshold 

} 

Block  { 

"74.5” 

BlockType 

Fen 

Name 

"Thrust  Tail  off' 

Position 

[195,  85,  255,  115] 

Expr 

} 

Block  { 

”u[2]*exp(-(u[l]- 

BlockType 

Outport 

Name 

"StOThrust” 

Position 

[410,165,430,  185] 

Port 

"1” 

OutputWhenDisabled 

held 

InitialOutput 

} 

Line  { 

”0” 

SrcBlock 

"Mux" 

SrcPort 

1 

DstBlock 

"Thrust  Tail  off' 

DstPort 

} 

Line  { 

1 

SrcBlock 

"Step" 

SrcPort 

1 

DstBlock 

"Mux" 

DstPort 

} 

Line  { 

2 

SrcBlock 

"Atmospere  Modell" 

SrcPort 

1 

DstBlock 

"StO  Thrust” 

DstPort 

1 

} 

Line  { 


SrcBlock 

"Altitude” 

SrcPort 

1 

DstBlock 

"Atmospere  Modell" 

DstPort 

1 

} 

Line  { 


SrcBlock 

"Switch” 

SrcPort 

1 

DstBlock 

"StOThrust" 

DstPort 

1 

} 

Line  { 


SrcBlock 

"Thrust  Tail  off’ 

SrcPort 

1 

DstBlock 

"Product” 

DstPort 

1 

} 


Line  { 


SrcBlock 

"StO  Thrust" 

SrcPort 

1 

Points 

[0,0] 

Branch  { 

Points 

[60,  0] 

DstBlock 

C-20 


DstPort 

} 

Branch  { 
Points 
DstBlock 
DstPort 

} 

} 

Line  { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 


3 


[0,  -95] 
2 


"Product" 


"Product" 

1 

[5,  0] 
"Switch" 

1 


Line  { 
SrcBlock 
SrcPort 
Points 
Branch  { 
DstBlock 
DstPort 

} 

Branch  { 
Points 
DstBlock 
DstPort 


"Time" 

1 

[15,  0] 

"Mux" 

1 


[0,  85] 
2 


"Switch" 


} 

} 

} 

} 

Block  { 

BlockType 

Name 

Ports 

Position 

NamePlacement 
ShowPortLabels 
System  { 

Name 

Location 

Open 

ScreenColor 
Block  { 
BlockType 
Name 
Position 
Port 

PortWidth 

SampIeTime 


Subsystem 
"Stl  Thrust" 

[2, 1,1,0, 0] 

[210,  172,  340,  203] 

alternate 

on 


"Stl  Thrust" 
[265,  309,  802,  476] 
off 
white 

Inport 
"Altitude" 
[110,  110,  130,  130] 
"1" 

"-1" 

"-1" 


} 

Block  { 
BlockType 
Name 
Position 
Port 

PortWidth 

SampIeTime 


Inport 

"Time" 

[85, 15, 105,  35] 
"2" 

"-1" 

"-l” 


} 

Block  { 
BlockType 
Name 
Ports 
Position 


EnablePort 

"Enable" 

[0,  0, 0,  0,  0] 
[35,15,55,  35] 


StatesWhenEnabling  held 
ShowOutputPort  off 


Block  { 


C-21 


BlockType 

Name 

Ports 

Position 

Inputs 

} 

Block  { 
BlockType 
Name 
Position 
Expr 

} 

Block  { 
BlockType 
Name 
Ports 
Position 
Inputs 

} 

Block  { 
BlockType 
Name 
Position 
Expr 

} 

Block  { 
BlockType 
Name 
Position 
Time 
Before 
After 

} 

Block  { 
BlockType 
Name 
Position 
Threshold 


Mux 

"Mux" 

[2,  1,0,  0,  0] 

[190,  17,  225,  53] 

"2" 


Fen 

"Pressure" 

[175,  105,235, 135] 

"101 325*exp((-l  .3728e-4)*u)" 


Product 

"Product" 

[2,  1, 0,  0,  0] 

[355,  27,  385,  58] 

"2" 


Fen 

"Stl  Thrust" 

[270,  105,  330, 135] 

"29803 2+(  1  -u/ 1 0 1 325)*  1 86854" 


Step 

"Step" 

[140,  30,  170,  60] 

-[147.7  154.7]" 

"[0  1]" 

"[1  01- 


Switch 

"Switch" 

[410,  65,  440,  95] 
"147.7" 


} 

Block  { 

BlockType 

Name 

Position 

Expr 

} 

Block  { 

BlockType 

Name 

Position 

Port 

OutputWhenDisabled 

InitialOutput 

} 

Line  { 

SrcBlock 
SrcPort 
Points 
Branch  { 

DstBlock 

DstPort 


Fen 

"Thrust  Tail  off 
[250,  20,310,  50] 

"u[2]*exp(-(u[l]-147.7)>" 


Outport 
"Stl  Thrust" 
[485, 70,  505,  90] 

T 

held 

"0" 


"Stl  Thrust" 

1 

[5,0] 

"Product" 

2 


} 

Branch  { 

Points  [55, 0] 

DstBlock  "Switch" 

DstPort  3 


} 


Line  { 


C-22 


SrcBlock 

SrcPort 

DstBlock 

DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 

) 

Line  { 
SrcBlock 
SrcPort 
Points 
Branch  { 
DstBlock 
DstPort 

} 

Branch  { 
Points 
DstBlock 
DstPort 

} 

} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 


"Switch" 

1 

"Stl  Thrust" 
1 


"Mux" 

1 

"Thrust  Tail  off 

1 


"Step" 

1 

"Mux" 

2 


"Thrust  Tail  off' 
1 

"Product" 

1 


"Product" 

1 

[5,  0] 

"Switch" 

1 


"Time" 

1 

[15,0] 

1 

[0,55] 

2 


"Mux" 


"Switch" 


"Pressure" 

1 

"Stl  Thrust" 
1 


} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

} 

} 

Block  { 

BlockType 

Name 

Position 

ShowName 

Threshold 

} 

Block  { 

BlockType 

Name 


"Altitude" 

1 

"Pressure" 

1 


Switch 

"Switch" 

[220,  25,  250,  55] 
off 

,,st3ig+69" 


Switch 

"Switch  1" 


Position 

[275,  60,  305,  90] 

ShowName 

off 

Threshold 

"st3ig" 

1 

Block  { 

BlockType 

Switch 

Name 

"Switch2" 

Position 

[330,  95,  360, 125] 

ShowName 

off 

Threshold 

i 

"80.7+st2ig" 

/ 

Block  { 

BlockType 

Switch 

Name 

"Switch3" 

Position 

[385,  130,415,160] 

ShowName 

off 

Threshold 

\ 

"st2ig" 

/ 

Block  { 

BlockType 

Switch 

Name 

"Switch4" 

Position 

[440,  165,  470,  195] 

ShowName 

off 

Threshold 

"154.7" 

} 

Block  { 


— l 

BlockType 

Switch 

Name 

"Switch5" 

Position 

[495,  200,  525,  230] 

ShowName 

off 

Threshold 

"81.5" 

} 

Block  { 


BlockType 

Switch 

Name 

,'Switch6" 

Position 

[140,  200,  170,  230] 

ShowName 

off 

Threshold 

"154.7" 

} 

Block  { 


BlockType 

Outport 

Name 

Position 

"Thrust" 

[550,  208,  580,  222] 

Port 

"1" 

OutputWhenDisabled 

held 

InitialOutput 

} 

Line  { 

"0" 

SrcBlock 

"Switch5" 

SrcPort 

1 

DstBlock 

"Thrust" 

DstPort 

} 

Line  { 

1 

SrcBlock 

"Switch" 

SrcPort 

1 

Points 

DstBlock 

[0,  25] 

"Switch  1" 

DstPort 

1 

} 

Line  { 


SrcBlock 

"Switch  1" 

SrcPort 

1 

Points 

[0,  25] 

DstBlock 

"Switch2" 

DstPort 

1 

} 

Line  { 

C-24 


Labels 

SrcBlock 

SrcPort 

Points 

DstBlock 

DstPort 

} 

Line  { 

SrcBlock 

SrcPort 

Points 

DstBlock 

DstPort 

} 

Line  { 

SrcBlock 

SrcPort 

Points 

DstBlock 

DstPort 

} 

Line  { 

Labels 
SrcBlock 
SrcPort 
Points 
Branch  { 

Points 
Branch  { 
DstBlock 
DstPort 

} 

Branch  { 
Points 
DstBlock 
DstPort 

} 

} 

Branch  { 

DstBlock 

DstPort 


[0,0] 

"Switch2" 

1 

[0,  25] 

"Switch3" 

1 


"Switch3" 

1 

[0,25] 

"Switch4" 

1 


"Switch4" 

1 

[0,  25] 

"Switch5" 

1 


[U0] 

"Motor  Off* 

1 

[5,0] 

[0,  -70] 

"Switch  1" 

3 


[0,  -55] 
"Switch" 

1 


"Switch3" 

3 


} 

i 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 


"St  2  Thrust" 

1 

"Switch2" 

3 


Line  { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 


"St  3  Thrust" 
1 

[0,  -20] 

"Switch" 

3 


"StO  Thrust" 

1 

[0,  -20] 

"Switch5" 

3 


} 

Line  { 

SrcBlock  "Time" 

SrcPort  1 

Points 
Branch  { 


[25,0] 


Points 
Branch  { 

DstBlock 
DstPort 

} 

Branch  { 

Points 
Branch  { 

DstBlock 
DstPort 

} 

Branch  { 

Points 
Branch  { 

DstBlock 
DstPort 

} 

Branch  { 

Points 
Branch  { 

DstBlock 
DstPort 

} 

Branch  { 

Points 
DstBlock 
DstPort 

} 

} 

} 

} 

} 

Branch  { 

Points 
Branch  { 

Points 
Branch  { 

Points 
Branch  { 

Points 
Branch  { 

Points 
Branch  { 

Points 
DstBlock 
DstPort 

} 

Branch  { 
DstBlock 
DstPort 

} 

} 

Branch  { 

DstBlock 
DstPort 

} 

} 

Branch  { 

DstBlock 
DstPort 

} 
i 

Branch  { 

DstBlock 
DstPort  2 

} 

) 

Branch  { 


[0,  30] 

"St  3  Thrust" 

1 


[0,  50] 

"St  2  Thrust" 

1 


[0,  95] 

"Switch6" 

2 


[0,  35;  80,  0] 

"StO  Thrust" 
2 


[0,  -55] 
"Stl  Thrust" 

2 


[100,  0] 

[0,  35;  55, 0] 

[0,  35;  55, 0] 

[0,  35;  55,  0] 

[0,  35;  55,  0] 

[0,  35] 

"SwitchS" 

2 


"Switch4" 

2 


"Switch3" 

2 


"Switch2" 

2 


"Switch  1" 


C-26 


DstBlock 

DstPort 


} 

} 

} 

Line  { 

SrcBlock 
SrcPort 
Points 
Branch  { 

DstBlock 

DstPort 


} 

Branch  { 

DstBlock 

DstPort 


} 

} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 


} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 


"Switch” 

2 


[100,0] 


"Switch6” 

1 


"Stl  Thrust" 

enable 


"StO  Thrust” 

enable 


"Constantl " 
1 

"Switch6" 

1 


"Constant" 

1 

"Switch6” 

3 


} 

Line  { 

SrcBlock 
SrcPort 
Points 
Branch  { 

DstBlock 

DstPort 

} 

Branch  { 

DstBlock 

DstPort 


"Altitude' 

1 

[135,  0] 

"Stl  Thrust" 

1 


"StO  Thrust" 

1 


} 

} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

} 


"Stl  Thrust" 
1 

"Switch4" 

3 


Block  { 
BlockType 
Name 
Ports 
Position 
ShowName 
Operator 

} 

Block  { 
BlockType 
Name 
Ports 
Position 
ShowName 
Operator 

} 

Block  { 


Trigonometry 
"TrigonometricVnFunction" 
[1,  1,0,  0,  0] 

[460,  315,  490,  345] 
off 
sin 


Trigonometry 

"Trigonometric\nFunctionl " 

[1,  1,0,  0,  0] 

[460,  360,  490,  390] 
off 
cos 


C-27 


BlockType 

Fen 

Name 

"local  gee" 

Position 

[275,  290,  335,  320] 

Orientation 

left 

Expr 

} 

Block  { 

"9.81  *(63781 35A2/(6378135+u)y 

BlockType 

Outport 

Name 

"Thrust" 

Position 

[695,  52,  725,  68] 

Port 

"1" 

OutputWhenDisabled 

held 

InitialOutput 

} 

Block  { 

"0" 

BlockType 

Outport 

Name 

"CG  Vel" 

Position 

[695,  163,725,  177] 

Port 

"2" 

OutputWhenDisabled 

held 

InitialOutput 

} 

Block  { 

"0" 

BlockType 

Outport 

Name 

"Density" 

Position 

[695,  17,  725,  33] 

Port 

"3" 

OutputWhenDisabled 

held 

InitialOutput 

\ 

"0" 

f 

Block  { 

BlockType 

Outport 

Name 

"Drag  Force" 

Position 

[320,  15,  340,  35] 

Port 

"4" 

OutputWhenDisabled 

held 

InitialOutput 

} 

Line  { 

"0" 

SrcBlock 

"Vibes" 

SrcPort 

1 

Points 

[230,  0] 

DstBlock 

"Muxl" 

DstPort 

} 

Line  { 

10 

SrcBlock 

"CG  Range" 

SrcPort 

1 

DstBlock 

,Terminator26" 

DstPort 

} 

Line  { 

1 

SrcBlock 

"Switch" 

SrcPort 

1 

Points 

[5,0] 

DstBlock 

"CG  Alpha" 

DstPort 

> 

Line  { 

1 

SrcBlock 

"Alpha  Dot" 

SrcPort 

1 

Points 

[0,  10] 

DstBlock 

"Switch" 

DstPort 

} 

Line  { 

1 

SrcBlock 

"Productl " 

SrcPort 

1 

C-28 


DstBlock 

DstPort 

} 

Line  { 
SrcBlock 
S  rcPort 
Points 
Branch  { 
DstBlock 
DstPort 

} 

Branch  { 
DstBlock 
DstPort 

} 

} 

Line  { 
SrcBlock 
SrcPort 
Points 
Branch  { 
Points 
DstBlock 
DstPort 

} 

Branch  { 
Points 
DstBlock 
DstPort 

} 

} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
Points 
Branch  { 
Points 
Branch  { 
DstBlock 
DstPort 

} 

Branch  { 
Labels 
Points 
DstBlock 
DstPort 

} 

} 

Branch  { 
Points 
DstBlock 
DstPort 

} 

} 

Line  { 


"CG  Range” 
1 


"Sum" 

1 

[0,  0] 

"TrigonometridnFunction" 

1 


"TrigonometricViFunctionl" 

1 


"local  gee" 

1 

[-10,0] 

[-115,0] 

"Mux” 

1 


[0,  -130] 

"Muxl" 

6 


"Mux" 

1 

"Alpha  Dot” 
1 


"Constant" 

1 

[30,  0] 
"Switch" 

3 


"Trigonometric\nFunction  1 " 
1 

[5,  0] 

[5,0] 

"Productl " 

1 


[2,  0] 

[0,  30;  -370,  0;  0,  -65] 
"Mux" 

4 


[0,  -105;  -230,  0;  0,  -75] 
"Muxl" 


8 


C-29 


SrcBlock 
SrcPort 
Points 
Branch  { 
DstBlock 
DstPort 

} 

Branch  { 
Points 
DstBlock 
DstPort 

} 

Branch  { 
Points 
DstBlock 
DstPort 

} 

} 

Line  { 
SrcBlock 
SrcPort 
Points 
Branch  { 
Points 
Branch  { 
DstBlock 
DstPort 

} 

Branch  { 
Points 
DstBlock 
DstPort 

} 


"Product" 

1 

[5,0] 

"CG  Altitude" 

1 


[0,  -60;  -300,  0;  0,  -60] 
"Muxl" 


9 


[0, 90;  -435,  0;  0,  -65] 
"Mux" 


5 


"Trigonometric\nFunction" 

1 

[0,  0] 

[10,0] 


"Product" 

2 


[0,  -50;  -240,  0;  0,  -95] 
"Muxl" 

7 


Branch  { 
Points 
DstBlock 
DstPort 

} 

} 

Line  { 
SrcBlock 
SrcPort 
Points 
Branch  { 
Points 
Branch  { 
DstBlock 
DstPort 


[0,  70;  -350,  0] 

"Mux" 

6 


"CG  Altitude" 

1 

[0,  -40] 

[-280,  0] 


"local  gee" 
1 


Branch  { 

Points  [-95, 0] 


Branch  { 

Points 
Branch  { 
Points 
DstBlock 
DstPort 

} 

Branch  { 
Points 
DstBlock 
DstPort 

} 

} 

Branch  { 

Points 

DstBlock 


[-125,  0] 

[0,  -100] 

"Thrust\nProgram" 

2 


[0,  35] 
"Mux" 
2 


[0,  -120] 

"Muxl" 


DstPort 

j 

} 

} 

Branch  { 

Points 

DstBlock 

DstPort 

} 

} 

Line  { 

SrcBlock 

SrcPort 

Points 

DstBlock 

DstPort 

} 

Line  { 

SrcBlock 

SrcPort 

DstBlock 

DstPort 

} 

Line  { 

SrcBlock 

SrcPort 

DstBlock 

DstPort 

} 

Line  { 

SrcBlock 
SrcPort 
Points 
Branch  { 

Points 

DstBlock 

DstPort 

} 

Branch  { 

Points 

DstBlock 

DstPort 

} 

} 

Line  { 

SrcBlock 

SrcPort 

DstBlock 

DstPort 


5 


[0,  -195] 

1 


"Atmospere  Model" 


"Pitch  Program" 
1 

[35,  0] 

"Sum" 

2 


"CG  Alpha" 
1 

"Sum" 

1 


"CG  Vel  Dot" 


"CG  Velocity" 
1 


"Drag" 

1 

H5,0] 
[0,  55] 
4 


"Muxl" 


[0,  -75] 


"Drag  Force" 


"Muxl" 

1 

"CG  Vel  Dot" 
1 


} 

Line  { 

SrcBlock 

SrcPort 

DstBlock 

DstPort 

} 

Line  { 

SrcBlock 
SrcPort 
Points 
Branch  { 

Points 
Branch  { 

Points 
Branch  { 

Points 
Branch  { 
Points 


"Mux2" 

1 

"Drag" 

1 


"CG  Velocity" 

1 

[30,  0] 

[0,0] 

[0,  145] 

[0,75] 

[0,  20;  -380,  0;  0,  -80] 


DstBlock 

DstPort 

} 

Branch  { 
DstBlock 
DstPort 

} 

} 

Branch  { 

DstBlock 

DstPort 


} 

} 

Branch  { 

Points 
Branch  { 

Points 

DstBlock 

DstPort 

} 

Branch  { 

DstBlock 

DstPort 

} 


} 

} 

Branch  { 
DstBlock 
DstPort 

} 


"Mux" 

3 


"Productl " 
2 


"Product" 

1 


[0,  -50;  -105,  0] 


"Muxl" 


[-120,  0] 

1 


"Mux2" 

2 


"CG  Vel" 


1 


Line  { 
SrcBlock 
SrcPort 
Points 
Branch  { 
Points 
DstBlock 
DstPort 

} 

Branch  { 
Points 
DstBlock 
DstPort 

} 

} 

Line  { 
SrcBlock 
SrcPort 
Points 
Branch  { 
DstBlock 
DstPort 

} 

Branch  { 
Points 
DstBlock 
DstPort 

} 

} 

Line  { 
SrcBlock 
SrcPort 
Points 
Branch  { 
DstBlock 
DstPort 

} 


"Thrust\nProgram" 

1 

[0,  -20;  15,0] 


[45,0] 

3 


"Muxl" 


[0,  -90] 

1 


"Thrust" 


"Atmospere  Model" 
1 

[-10,  0] 

"Mux2" 

1 


[0, -65] 
1 


"Density" 


"Time" 

1 

[70,  0] 

"Thrust\nProgram" 

1 


Branch  { 

Points 

Branch  { 
Points 
DstBlock 
DstPort 

} 

Branch  { 
Points 
DstBlock 
DstPort 

} 

} 

} 

Line  { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 

} 


[0,  210] 

[150,  0;  0,  -10] 

"Switch" 

2 


[0,  80] 

"Pitch  Program" 
1 


"Sys  Mass" 

1 

[180,  0;  0, 15] 
"Muxl" 

2 


} 

Block  { 

BlockType 

Name 

Ports 

Position 

ShowPortLabels 
System  { 

Name 

Location 

Open 

ScreenColor 
Block  { 
BlockType 
Name 
Position 
Port 

PortWidth 

SampleTime 

} 

Block  { 
BlockType 
Name 
Position 
Port 

PortWidth 

SampleTime 


Subsystem 

"Mass  Matrix  \nlnverter" 
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[880,  425,  935,  445] 
"massl" 


Step 

"Neg  Nozzle” 

[400,  455,  430,  485] 

*'[0;81  -5;st2ig;st3ig] ,f 
”[-795.3;-235.4;-l  18.7;-52.8],r 
”[0;0;0;0]" 


Constant 
"Payload  Mass” 
[355,  50,  410,  70] 

"pirn” 


Step 

"Pos  Nozzle” 

[400,415, 430,  445] 
alternate 

"  [8 1 .5 ;  1 54.7  ;st2ig+80.7;st3ig+69]  ” 
”[795.3;235.4;118.7;52.8]" 

”[0;0;0;0]" 


Block  { 
BlockType 
Name 
Position 
Time 
Before 
After 


Step 

”St0  Mass" 
[425, 65,  455,  95] 
”81.5” 

”0” 

"-42 11" 


} 

Block  { 
BlockType 
Name 
Position 
Time 
Before 
After 

} 

Block  { 
BlockType 
Name 
Position 
Time 
Before 
After 

} 

Block  { 
BlockType 
Name 
Position 
Time 
Before 
After 

} 

Block  { 
BlockType 
Name 
Position 
Time 
Before 
After 

} 

Block  { 
BlockType 


Step 

"Str 

[705,  330,  735,  360] 
"81.5" 

"14020” 

"0" 


Step 

"Stl  Mass" 
[475,  85,  505,  115] 
"154.7" 

"0" 

”-1868" 


Step 

"Str 

[665,  360,  695,  390] 
"154.7" 

”3370” 

"0” 


Step 

"St2  Mass" 
[525, 105,  555,  135] 
"st3ig" 

”0  ” 

”-345” 


Step 
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Name 

Position 

Time 

Before 

After 

} 

Block  { 
BlockType 
Name 
Position 
Time 
Before 
After 

} 

Block  { 
BlockType 
Name 
Ports 
Position 
Inputs 

} 

Block  { 
BlockType 
Name 
Ports 
Position 
Orientation 
Inputs 

} 

Block  { 
BlockType 
Name 
Position 
ShowName 
Threshold 

} 

Block  { 
BlockType 
Name 
Position 
ShowName 
Threshold 

} 

Block  { 
BlockType 
Name 
Position 
ShowName 
Threshold 

} 

Block  { 
BlockType 
Name 
Position 
ShowName 
Threshold 

} 

Block  { 
BlockType 
Name 
Position 
ShowName 
Threshold 

} 

Block  { 
BlockType 
Name 
Position 


"St3" 

[700,  390,  730,  420] 
"st3ig" 

"985" 

"100" 


Step 

"St3  Mass" 
[620,  150,  650, 180] 

"st3ig+129" 

"0" 

"-203" 


Sum 

"Sum" 

[8,  1,0,  0,0] 
[530, 406,  565,  494] 
"++++++++" 


Sum 

"Suml" 

[2,  1, 0,  0,  0] 
[365,  342,  395,  373] 
left 


Switch 
"  Switch  1" 
[375,  85, 405,115] 
off 

"st3ig+69" 


Switch 
"Switch2" 
[425,  110,  455,  140] 
off 

"st3ig" 


Switch 
"Switch3" 
[475,  135,  505, 165] 
off 

"st2ig+80.7" 


Switch 
"Switch4" 
[625,  210,  655,  240] 
off 

"81.5" 


Switch 
"Switch5" 
[525, 160,  555,  190] 
off 

"st2ig" 


Switch 
"Switch6" 
[575,  185, 605,215] 


ShowName  off 

Threshold  "154.7" 


} 

Block  { 
BlockType 
Name 
Ports 
Position 
Inputs 

} 

Block  { 
BlockType 
Name 
Position 
Value 


Sum 

"System  Mass" 
[7, 1,0,  0,  0] 
[760, 127,  790,  233] 
"+++++++" 


Constant 

"mdot’s" 

[30,  144,  270,  166] 

"[-11.3333;-37 .48451  ;~166.0109;-598.8834;0]" 


} 

Block  { 

BlockType 

Name 

Position 

Port 

OutputWhenDisabled 

InitialOutput 

} 

Block  { 

BlockType 

Name 

Position 

Orientation 

NamePlacement 

Port 

OutputWhenDisabled 

InitialOutput 


Outport 
"Tot  Mass" 
[970, 170,  990,  190] 
"1" 
held 
"0" 


Outport 
"Lower" 
[320,  350,  340,  370] 
left 

alternate 
"2 " 
held 
"0" 


Block  { 

BlockType 

Name 

Position 

Orientation 

Port 

OutputWhenDisabled 

InitialOutput 

} 

Block  { 

BlockType 

Name 

Position 

Port 

OutputWhenDisabled 

InitialOutput 

} 

Block  { 

BlockType 

Name 

Position 

Port 

OutputWhenDisabled 

InitialOutput 

} 

Block  { 

BlockType 

Name 

Position 

Port 

OutputWhenDisabled 

InitialOutput 

} 

Line  { 


Outport 
"Upper" 
[320,  295,  340,  315] 
left 
"3" 
held 
"0" 


Outport 

"Inert" 

[960,380,  980,  400] 
"4" 
held 
"0" 


Outport 
"Payload" 
[955, 270,  975,  290] 
"5" 
held 
"0" 


Outport 

"SCI" 

[960,  425,  980,  445] 
"6" 
held 
"0" 
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SrcBlock 

SrcPort 

DstBlock 

DstPort 


"Fcnl" 

1 

"Payload’ 

1 


} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 


"Fairing" 

1 

"Inert\nMass" 

4 


} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 


"St2" 

1 

"Inert\nMass" 

2 


} 

Line  { 

SrcBlock 
SrcPort 
Points 
Branch  { 

Points 
Branch  { 

Points 
Branch  { 

Points 
Branch  { 
Points 
Branch  { 
Points 
DstBlock 
DstPort 


"Time" 

1 

[190,  0] 

[0,  25;  50,  0] 
[0,  25;  50,  0] 


[0,  25;  50, 0] 
[0,  25;  50,  0] 


[0,  25] 
2 


"Switch4" 


} 

Branch  { 

DstBlock  "Switch6" 

DstPort  2 


} 

} 

Branch  { 
DstBlock 
DstPort 

} 

} 

Branch  { 

DstBlock 

DstPort 

} 

} 

Branch  { 

DstBlock 

DstPort 


} 

Branch  { 

DstBlock 

DstPort  2 


"Switch5" 

2 


"Switch3" 

2 


"Switch2" 

2 


"Switchl" 


} 

} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 


"Switchl1 

1 

"Switch2’ 

1 


"Switch2' 

1 


DstBlock 

DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 


"Switch3" 

1 


"Switch3" 

1 

"SwitchS" 

1 


"SwitchS" 

1 

"Switch6" 

1 


"Switch6" 

1 

"Switch4" 

1 


"Switch4" 

1 

"Mass  Flow" 

1 


"Mass  Flow" 

1 

"System  Mass" 
7 


"St3  Mass" 

1 

[5,  0;  0,  45] 
"System  Mass" 
6 


Line  { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 

i 

Line  { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 


"Fairing\nMass" 

1 

[55,  0;  0,  55] 
"System  Mass" 

5 


"St2  Mass" 

1 

[110,  0;  0,  60] 
"System  Mass" 
4 


"Stl  Mass" 

1 

[165, 0;  0, 65] 
"System  Mass" 
3 


"StO  Mass" 

1 

[220,  0;  0,  70] 
"System  Mass" 
2 


Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 


"mdot's" 

1 

’’Demux" 

1 


"Demux" 

1 

[20, 0;  0,-15] 
"Switch  1" 

3 


"Demux" 

2 

[20,  0;  0,  20] 

"Switch3" 

3 


"Demux" 

3 

[15, 0;  0, 55] 
"Switch6" 

3 


} 

Line  { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 

} 

Line  { 

Labels 
SrcBlock 
SrcPort 
Points 
Branch  { 
Points 
Branch  { 
Points 
DstBlock 
DstPort 

} 

Branch  { 
DstBlock 
DstPort 

} 

} 

Branch  { 
DstBlock 
DstPort 

} 

} 

Line  { 
SrcBlock 
SrcPort 
Points 
Branch  { 
DstBlock 
DstPort 

} 

Branch  { 
Points 
DstBlock 
DstPort 


"Demux" 

4 

[10,  0;  0,  65] 
"Switch4" 

3 


[1,0] 

"Demux" 

5 

[25,  0] 

[0,  -50] 

[0,  -45] 

"Switchl" 

1 


"Switch2" 

3 


"Switch5" 

3 


"System  Mass" 

1 

[25,0] 

"Tot  Mass" 

1 


[0,  100] 

1 


"Active  Stage  Mass" 


} 

} 

Line  { 
SrcBlock 
SrcPort 
Points 
Branch  { 
Points 
DstBlock 
DstPort 

} 

Branch  { 
Points 
Branch  { 
DstBlock 
DstPort 

} 

Branch  { 
Points 
DstBlock 
DstPort 

} 


} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 


"Payload  Mass" 

1 

[270,  0] 


[0,75] 

1 


"System  Mass” 


[150,  0;  0,  220] 


"FcnlM 


1 


[0,  15] 

"Active  Stage  Mass" 
2 


"Pos  Nozzle" 
1 

"Demux  1" 

1 


"Neg  Nozzle" 

1 

"Demux2" 

1 


"Demux  1" 
1 

"Sum" 

1 


"Demux  1" 
2 

"Sum" 

2 


"Demux  1" 
3 

"Sum" 

3 


"Demux  1" 
4 

"Sum" 

4 


"Demux2" 

1 

"Sum" 

5 


"Demux2" 
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SrcPort 

DstBlock 

DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 


2 

"Sum" 

6 


"Demux2" 

3 

"Sum" 

7 


} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 


"Demux2" 

4 

"Sum" 

8 


"Active  Stage  Mass" 
1 

"Fen" 

1 


} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
Points 
Branch  { 
DstBlock 
DstPort 

} 

Branch  { 
Points 
DstBlock 
DstPort 


"Suml" 

1 

"Lower" 

1 


"St3" 

1 

"Inert\nMass" 

3 


"Stl" 

1 

"Inert\nMass" 

1 


"Massl" 

1 

"SCI" 

1 


"Sum" 

1 

[25,  0;  0,  -85] 

"Suml" 

2 


[0,  -40] 
4 


"Active  Stage  Mass" 


} 

} 

Line  { 

SrcBlock  "Inert\nMass" 

SrcPort  1 

Points  [15,0] 

Branch  { 

DstBlock  "Inert" 

DstPort  1 

} 
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Branch  { 
Points 
DstBlock 
DstPort 

} 

} 

Line  { 
SrcBlock 
SrcPort 
Points 
Branch  { 
DstBlock 
DstPort 

} 

Branch  { 
DstBlock 
DstPort 

} 

.  } 

} 

} 

Block  { 

BlockType 

Name 

Ports 

Position 

ShowPortLabels 

MaskType 

MaskDescription 

MaskPromptString 


MaskStyleString 
Mask  Variables 

MasklconFrame 

MasklconOpaque 

MasklconRotate 

MasklconUnits 

MaskValueString 

System  { 

Name 

Location 

Open 

ScreenColor 
Block  { 
BlockType 
Name 
Position 
Port 

PortWidth 

SampleTime 

} 

Block  { 
BlockType 
Name 
Position 
Port 

PortWidth 

SampleTime 

} 

Block  { 
BlockType 
Name 
Position 
Port 

PortWidth 

SampleTime 


[0,  -80] 

"Active  Stage  Mass" 
3 


"Fen" 

1 

[-25,0] 

"Upper" 

1 


"Suml" 

1 


Subsystem 
"Noise\n  Generator" 

[3,  1,0,0,  0] 

[440,  58,  560,  112] 
on 

"Noise  Generator" 

"Set  band  limited  noise  power." 

"Aero  Noise  PowerIResonant  Bum  FrequencyIResonan" 
"t  Bum  AmplitudelCombustion  Noise  PowerlNoise  " 
"Filter  NumeratorINoise  Filter  Denominator" 

"edit, edit, edit, edit, edit, edit" 

"  wind=  @  1  ;freq=  @  2  ;respo  wer=  @  3  ;power=  @  4  ;bp=@  5  ;ap=  @  6 " 

on 

on 

none 

autoscale 

"5.0e5l50l5e4!le2lbplap" 

"Noise\n  Generator" 

[371,513,  932,883] 
off 
white 

Inport 

"Thrust" 

[20,215,  40,  235] 

(i 

"-r 

"-i" 


Inport 
"Velocity" 
[20,  250,  40,  270] 

"2" 

M" 

M" 


Inport 
"Density" 
[20,  285,  40,  305] 

"3" 

"-l" 

"  - 1 " 
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} 

Block  { 
BlockType 
Name 
Ports 
Position 
SourceBlock 
SourceType 
Cov 
Ts 
seed 


Reference 

"Band-Limited\nWhite  Noise" 

[0,  1, 0,  0,  0] 

[20, 150,  50,  180] 

"simulink/Sources/Band-Limited\n  White  Noise" 
"Continuous  White  Noise." 

"[1  wind]" 

"1/400" 

"[12312  34985]" 


Block  { 
BlockType 
Name 
Ports 
Position 
NamePlacement 
ShowPortLabels 
System  { 

Name 

Location 

Open 

ScreenColor 
Block  { 
BlockType 
Name 
Position 
Port 

PortWidth 

SampleTime 

} 

Block  { 
BlockType 
Name 
Position 
Port 

PortWidth 

SampleTime 

} 


Subsystem 
"Combustion  Noise" 
[2,  1,  0,  0,  0] 

[210, 132,  290, 183] 
alternate 


on 

"Combustion  Noise" 

[576,  28,  871,  203] 
off 


white 

Inport 

"Noise  In" 


"1" 


[20,  55,  40,  75] 
"-l" 


"-1" 


"Thrust" 

"2" 

"-1" 


Inport 

[20,  115,40,  135] 
"-l" 


Block  { 

BlockType 

Name 

Position 

Expr 

} 

Block  { 

BlockType 

Name 

Ports 

Position 

Inputs 

} 

Block  { 

BlockType 

Name 

Position 

Port 

OutputWhenDisabled 

InitialOutput 

} 

Line  { 

SrcBlock 

SrcPort 

Points 

DstBIock 

DstPort 

} 


Fen 

"Fen" 

[160,  60,  220,  90] 
"u[l]*u[2]/1650000" 


Mux 

"Mux" 

[2,  1,0,0,  0] 

[105,  57, 140,  93] 

"2" 


Outport 
"Noise  Out" 

[240,  68,  270,  82] 
"l" 
held 

"0" 


"Thrust" 

1 

[45,0] 

"Mux" 

2 
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Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line  { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

} 

} 

Block  { 
BlockType 
Name 
Ports 
Position 
Outputs 

} 

Block  { 
BlockType 
Name 
Position 
Gain 


"Noise  In" 

1 

"Mux" 

1 


"Mux" 

1 

"Fen" 

1 


"Fen" 

1 

"Noise  Out" 

1 


Demux 
"Demux" 
[1,2,  0,  0,  0] 
[80,  147,  115,  183] 

"2" 


Gain 

"Gain" 

[320,  144,380,  176] 
"power" 


Block  { 
BlockType 
Name 
Ports 
Position 
Inputs 

} 

Block  { 
BlockType 
Name 
Ports 
Position 
Inputs 

} 

Block  { 
BlockType 
Name 
Ports 
Position 
ShowName 
Inputs 

} 

Block  { 
BlockType 
Name 
Position 
Amplitude 
Frequency 
Phase 

SampleTime 

} 

Block  { 
BlockType 
Name 
Position 


Mux 

"Mux" 

[2,  1, 0,  0,  0] 

[455,  132,  490, 168] 

”2" 


Product 

"Product" 

[2, 1,0,  0,0] 

[185,32,215,  63] 

"2" 


Product 
"Product  1" 
[2,  1,  0, 0,  0] 
[195,  257,  225,  288] 
off 

"2" 


Sin 

"Resonant  Bum" 
[125,  60, 155,  90] 

"respower" 

"(freq)*2*pi" 

"0" 

"0" 


Step 
"Step" 
[125,15,155,  45] 


’[81.5] 


Time 

Before 

After 

} 

Block  { 
BlockType 
Name 
Ports 
Position 
Inputs 

} 

Block  { 
BlockType 
Name 
Position 
Numerator 
Denominator 

} 

Block  { 
BlockType 
Name 
Position 

NamePlacement 

Numerator 

Denominator 

} 

Block  { 
BlockType 
Name 
Ports 
Position 
ShowPortLabels 
System  { 

Name 

Location 

Open 

ScreenColor 

Block  { 
BlockType 
Name 
Position 
Port 

PortWidth 

SampleTime 

} 

Block  { 
BlockType 
Name 
Position 
Port 

PortWidth 

SampleTime 

} 

Block  { 
BlockType 
Name 
Position 
Expr 

} 

Block  { 
BlockType 
Name 
Position 
Gain 

} 

Block  { 
BlockType 
Name 


"[!]" 


[0]" 


Sum 

"Sum" 

[2,  1, 0,  0,  0] 

[400,71,430,  104] 

"++" 


TransferFcn 
"Transfer  Fen" 
[140,  165,  175,  205] 

"bp" 

"ap" 


TransferFcn 
"Transfer  Fen  1" 
[140, 125, 175, 165] 
alternate 
"bp" 

"ap" 


Subsystem 

"Transonic  Scale  Factor" 

[2,  1, 0,  0,  0] 

[75,  251, 175,  304] 


"Transonic  Scale  Factor" 

[125,  273,  644,  597] 
off 

white 


"Velocity" 

"1" 

"-l" 


Inport 

[15,  88,  45,  102] 
"-l" 


"Density" 

"2" 

"-r 


Inport 

[15,  223,  45,  237] 
"-1" 


Fen 

"Dynamic  Pressure" 

[150,205,210,235] 

"(u[2]*u[l]A2)/2" 


"Gain" 

"12" 


Gain 

[245,  136,  285,  164] 


Fen 

"Mach  (V/c)  " 


C-49 


Position 

Expr 

} 

Block  { 
BlockType 
Name 
Position 
Expr 

} 

Block  { 
BlockType 
Name 
Ports 
Position 
Inputs 

} 

Block  { 
BlockType 
Name 
Ports 
Position 
Inputs 


[80,  80,  140,  110] 
"u/(L4*287.05*249)A0.5" 


Fen 

"Mach  Effects" 

[250,  90,310,  120] 
"l/(0.5+(abs(l-u[l]A2))A0.5)" 


Mux 

"Mux" 

[2,  I,  0,  0,  0] 

[185,  87,220,  123] 

"2" 


Mux 

"Muxl" 

[2,1,0,  0,  0] 

[95,  202,130,  238] 

”2" 


Block  { 
BlockType 
Name 
Ports 
Position 
Inputs 

} 

Block  { 
BlockType 
Name 
Position 
Expr 

} 

Block  { 
BlockType 
Name 
Ports 
Position 
Inputs 

} 


Product 

"Product" 

[2,  1,  0,  0,  0] 

[415,  112, 445, 143] 

"2" 


Fen 

"Scale  Pressure" 

[75,  135,  135,  165] 

"u/101325" 


Sum 

"Sum" 

[2,1,0,  0,  0] 

[345,  102,  375, 133] 

"++" 


Block  { 

BlockType 

Name 

Position 

Port 

OutputWhenDisabled 

InitialOutput 

} 

Line  { 

SrcBlock 

SrcPort 

Points 

DstBlock 

DstPort 


"Scale" 

"1" 

held 

"0" 


Outport 

[480,  123,510,  137] 


"Dynamic  Pressure" 
1 

[0,  -35; -155,0] 

"Scale  Pressure" 

1 


} 

Line  { 

SrcBlock 
SrcPort 
Points 
Branch  { 

Points 

DstBlock 

DstPort 

} 


[5,0] 


"Velocity" 

1 


"Muxl" 


[0,115] 

1 


Branch  { 


DstBlock 

DstPort 

} 

} 

Line  { 

SrcBlock 

SrcPort 

DstBlock 

DstPort 

} 

Line  { 

SrcBlock 

SrcPort 

DstBlock 

DstPort 

} 

Line  { 

SrcBlock 

SrcPort 

DstBlock 

DstPort 

} 

Line  { 

SrcBlock 
SrcPort 
Points 
Branch  { 

Points 

DstBlock 

DstPort 

} 

Branch  { 

DstBlock 

DstPort 

} 

} 

Line  { 

SrcBlock 

SrcPort 

DstBlock 

DstPort 

} 

Line  { 

SrcBlock 

SrcPort 

Points 

DstBlock 

DstPort 

} 

Line  { 

SrcBlock 

SrcPort 

Points 

DstBlock 

DstPort 

} 

Line  { 

SrcBlock 

SrcPort 

DstBlock 

DstPort 


"Mach  (V/c) " 
1 


"Muxl" 
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